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APPENDIX  I 

Detailed  Workplan 

Predicting  the  Behavior  of  Orimulsion 
Spilled  on  Water 


EXECUTIVE  SUMMARY 


Orimulsion,  an  alternative  fuel  for  power  generation,  is  currently  being  used  in  countries  such  as  Canada, 
Italy,  Denmark,  Japan  and  China  and  has  been  tested  in  the  U.S.  in  Florida  and  Illinois.  To  each  of  these 
locations,  Orimulsion  is  being  shipped  in  ocean  going  tankers.  In  order  to  assess  the  required 
countermeasures  in  case  of  a  spill,  a  thorough  understanding  of  its  behavior  when  released  in  water  is 
required."  Bench  top  and  flume  tank  studies  to  date  have  demonstrated  the  complexities  surrounding  its 
behavior  in  the  event  of  a  release.  For  example,  previous  studies  have  shown  that  receiving  water  salinity 
will  markedly  effect  Orimulsion’s  behavior.  However,  the  impact  and  interactions  of  potentially 
important  independent  variables  other  than  salinity,  e.g.,  water  temperature,  Orimulsion  concentration, 
particulate  load,  particulate  type  and  turbulence  of  the  water  are  poorly  known.  As  a  result,  the  need  exists 
to  rigorously  evsduate  the  behavior  of  Orimulsion  under  a  wide  range  of  environmental  conditions.  This 
wori^lan  describes  a  joint  study  devised  by  Battelle  Memorial  Institute,  Bitor  America  Corp.  Bitumenes 
Orinoco,  S.A.,  Intevep  S.A.,  United  States  Coast  Guard,  and  National  Oceanogr^hic  and  Atmospheric 
Administration  and  co-fimded  by  Bitor  America  Corp.  and  the  United  States  Coast  Guard  Research  and 
Development  Center. 

The  workplan  consists  of  two  principal  phases  -  a  benchtop  study  phase  and  a  meso-scale  flume  tank 
study  phase.  The  benchtop  phase  utilizes  design  of  experiment  (DOE)  principles  in  order  to  reduce  the 
number  of  experiments  necessary  to  achieve  a  statistically  valid  result.  These  experiments  will  be 
conducted  in  4L  glass  vessels.  Initially,  the  benchtop  experiments  will  include  a  12-run  Plackett-Burman 
screening  design  phase  in  which  the  gross  behavior  of  Orimulsion  will  be  evaluated  rmder  a  wide  range  of 
environmental  conditions  represented  by  the  six  independent  variables  of  interest  (salinity,  temperature, 
Orimulsion  concentration,  particulate  loading,  particulate  type,  and  turbulence).  The  screening  design 
will  determine  if  any  of  &e  six  independent  variables  has  no  impact  on  Orimulsion’s  behavior,  as 
measured  by  particle  growth  and  its  ‘compartmentalization’  (i.e.,  does  it  float,  sink,  remain  suspended, 
dissolve,  or  stick  to  walls)  after  24  hours.  Four  high  Orimulsion  concentration  (50,000  mg/L) 
experiments  will  also  be  performed  during  this  phase.  Those  variables  that  influence  Orimulsion’s 
behavior  will  subsequently  be  evaluated  and  quantified  in  a  3-level  fectorial,  fece-centered  cube  design 
phase  that  will  include  up  to  47  additional  experiments.  The  objectives  of  the  3-level  design  are  to  (1) 
quantify  the  influence  of  the  individual  independent  variables  and  interaction  between  independent 
variables  on  Orimulsion’s  behavior  and  (2)  help  focus  the  conditions  under  which  the  meso-scale  flume 
tank  studies  will  be  conducted  in  the  second  phase  of  the  study. 

In  the  meso-scale  flume  tank  study  phase,  a  limited  number  of  experiments  will  be  conducted  in  Battelle’s 
4000L  flume  tank  under  enviromnental  conditions  identified  as  critical  in  the  benchtop  phase.  The  results 
of  the  flume  studies  will  be  used  to  groimd  tmth  the  benchtop  studies  under  conditions  that  more  closely 
mimic  open  water  conditions.  Prior  to  the  Orimulsion  release  experiments,  an  attempt  will  be  made  to 
quantify  the  turbulent  diffusive  energy  within  the  flume  tank  under  varying  wind  speed,  wave  height,  and 
wave  fiequency  conditions.  This  information  is  necessary  in  order  to  project  the  results  fi-om  the  flume 
experiments  to  open  water  conditions.  Realizing  the  shortcomings  of  any  laboratory  experiments  in 
mimicking  nature,  the  ultimate  objectives  of  this  study  are  to  improve  the  understanding  of  Orimulsion’s 
behavior  when  spilled  on  water  and  to  provide  practical  guidelines  for  responders  in  the  event  of  an 
Orimulsion  spill  under  a  variety  of  conditions. 
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1.  INTRODUCTION 

In  the  event  of  an  oil  spill  at  sea,  most  oils  will  float  and  rapidly  begin  to  undergo  feiriy  predictable 
modifications  to  their  physical  and  chemical  properties  collectively  known  as  weathering.  Most  oil  spill 
contingency  measures  currently  available  have  been  developed  for  floating  oils  and  their  respective 
emulsions.  Recently,  however,  there  is  increased  interest  in  predicting  the  behavior  of  a  relatively  new, 
petroleum-derived  fuel  known  as  Orimulsion. 

Orimulsion  is  a  fuel  that  is  an  emulsion  made  of  ^proximately  73%  natural  bitumen  and  27%  ficsh  water 
that  is  of  great  interest  to  electric  power  utilities  because  of  its  competitive  cost  and  pricing  structure. 
Orimulsion  is  in  commercial  use  in  other  coimtries  (e.g.,  Denmark,  Canada,  Japan,  Italy,  China, 
Lithuania,  Barbados)  but  is  not  in  use  in  the  United  States.  At  tire  present  time,  Florida  Power  and  Light 
Co.  has  proposed  to  transport  Orimulsion  through  Tampa  Bay  for  their  Manatee  plant.  The  permit  for 
that  operation  is  being  reviewed  by  the  Governor  of  Florida  and  his  Cabinet.  If  ^proved,  FPL  expects  to 
bring  Orimulsion  into  Tampa  Bay  in  two  years.  In  addition  to  FPL  there  are  thought  to  be  severd  other 
East  Coast  oil-fired  power  plants  considering  conversion  to  Orimulsion.  As  the  potential  for  transport  of 
Orimulsion  in  U.S.  water  is  eminent,  the  need  exists  to  better  understand  tire  fete  of  Orimulsion  if  it  were 
to  be  accidentally  spilled. 


2.  SCOPE  OF  STUDY 

Orimulsion  is  classified  as  a  Group  V  product  along  with  ofeer  heavy  oils  called  low  API  oils  (LAPIO). 
Various  benchtop  scale  studies  of  Orimulsion  and  other  LAPIO  fuels,  in  particular  fuel  oil  #6,  have  been 
performed  (Brown  et  al.,  1995;  Febres  et  al.,  1995;  Potter,  1995).  Battelle’s  1997  meso-scale  (i.e.,  flume 
tank)  study  of  LAPIO  fuels  included  Orimulsion  (Ostazeski  et  al.,  1997a).  During  this  study  14.2  L  of 
Orimulsion  was  released  into  a  4000  L  flume  tank  (3542  mg  Orimulsion/L  water)  containing,  on  separate 
occasions,  freshwater  (0  °/oo)  and  salt  water  (32  7oo)-  It  was  determined  that  the  behavior  of  Orimulsion 
100  when  spilled  in  water  is  very  different  fi’om  other  Group  V  products  and  as  a  result,  it  presents  a 
variety  of  very  challenging  response  scenarios. 

The  receiving  water  salinity  was  determined  to  be  a  critical  factor  effecting  Orimulsion’s  behavior.  In  the 
flesh  water  flume  test,  the  Orimulsion  was  observed  to  disperse  in  a  black  cloud  of  fine  droplets  that 
remained  in  suspension  throughout  the  120-hour  experiment.  In  the  sea  water  flume  test,  Orimulsion 
initially  dispersed  throughout  the  water  colurtm  but  within  15  minutes,  a  heavy  viscous  product  composed 
of  coalesced  bitumen  particles  began  to  float  on  the  surface.  The  surface  slick  continued  to  ‘grow’  by 
coalescence  throughout  the  120-hour  experiment.  The  complex  interaction  of  bitumen  particle  buoyancy, 
turbulence,  and  reduced  effectiveness  and  concentration  of  the  polyethoxylated  nonylphenol  surfectant  in 
sea  water  was  believed  to  allow  the  bitumen  particles  to  either  (1)  coalesce  and  rise  to  the  surface  forming 
the  viscous  surface  slick,  (2)  rise  to  the  surfece  and  then  coalesce,  or  (3)  a  combination  of  these  processes. 

At  this  point,  the  behavior  of  Orimulsion  in  salinities  representative  of  brackish  estuaries  is  uncertain. 
Similarly,  the  effect  that  other  potentially  important  variables,  e.g.,  water  temperature,  Orimulsion 
concentration,  particulate  load,  particulate  type,  and  turbuJpnce  of  the  water  (i.e.,  energetics),  have  on 
Orimulsion’s  behavior  are  not  well  known.  Therefore  in  the  event  of  a  release,  Orimulsion  could  be 
found  to  predominantly  float,  or  sink,  or  remain  suspended  in  tire  water  column.  The  need  to  predict  the 
fete  of  Orimulsion  imder  a  variety  of  environmental  conditions  is  of  critical  importance  for  containment, 
recovery  and  spill  modeling. 

This  workplan  outlines  Battelle's  proposed  approach  to  gaining  an  understanding  of  the  important 
parameters  that  affect  the  behavior  of  Orimulsion  when  spilled  on  water.  The  objective  of  the  study  is  to 
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investigate  the  weathering  characteristics  of  Orimulsion  in  water,  particulariy  as  it  relates  to  the  ultimate 
fate  of  spilled  Orimulsion.  The  project  was  proposed  to  consist  of  three  phases,  namely; 

□  Phase  I  -  Design  of  Behavior  Characterization  Experiments 

□  Phase  n-  Evaluation  ofFactors  Effecting  Orimulsion’s  Fate  in  Water 

•  Phase  DA  -  design  and  conduct  an  ^piopriate  number  of  benchtop 
e;q)eriments  to  evaluate  the  mnnerous  independent  variables  effecting 
Oiimulsion’s  fete  based  on  design  of  ejqperiment  (DOE)  principles 

•  Phase  DB  -  evaluate  bench-top  data  using  DOE  principles  to  identify  the 
most  significant  fectors  effecting  the  fete  of  Orimulsion. 

•  Phase  nC  -  further  investigate  and  quantify  flume  tank  energetics  so  that 
Phase  in  can  be  conducted  under  more  representative  wave  and  wind 
conditions  (than  previous  flume  experiments). 

□  Phase  in  -  Perform  meso-scale  flume  ejqperiments  that  focus  on  conditions 
which  have  the  greatest  effect  on  Orimulsion’s  fete  in  water  (as  based  on  the 
results  of  Phase  11)  and  prepare  a  detailed  comprehensive  report  of  the  entire 
study. 

This  report  comprises  the  deliverable  for  Phase  I  of  this  study.  The  objective  of  this  phase  is  to  propose  a 
detailed  course  of  study  for  Phases  n  and  in.  Upon  review  and  concurrence  of  this  course  of  study  by 
both  the  U.S.  Coast  Guard  (USCG)  and  BITOR  America  Corp.  (BITOR),  Phase  n  will  be  initiated. 


3.  PHASE  II  -  EVALUATION  OF  FACTORS  EFFECTING  THE  FATE  OF 
ORIMULSION 

3.1  Introduction  to  Experimental  Design 

The  independent  variables  identified  as  potentially  impading  fee  fete  of  Orimulsion  spilled  in  water  were  discussed  and 
agreed  iqxm  at  a  meeting  held  at  our  Duxbury  fecility  on  Oct.  31,  1997  which  was  attended  by  BITOR, 
NOAA,  and  USCG  representatives.  These  variables  and  their  ranges  of  interest  included; 


Independent  Variables  Ranee 

0-35  ppt 


water  salinity 
water  temperature 
Orimulsion  cone, 
particulate  loading 
particulate  type 
water  energetics 


0-30  X 

1000-17000  mg/L 
none  to  heavy 
calcite/kaolinite 
calm  to  heavy  seas 


Selected  Values  of  Interest 

0,  17,  35  ppt 
3,  15,  25  X 

1000-16000  mg/L;  50,000  mg/L" 
none,  20-120  mg/L 
calcite/kaolirrite 
static  to  dynamic 


With  six  variables  and  19  degrees  of  fieedom  an  ad  hoc  approach  (i.e.,  changing  only  one  variable  at  a 
time  while  holding  the  others  constant)  would  require  648  separate  experiments.  The  costs  and 
interpretation  of  such  an  approach  would  be  impractical.  Furthermore,  tlris  ‘one  variable  at  a  time 
^proach’  ignores  the  interaction(s)  that  may  exist  among  the  individual  variables.  Therefore,  an 
alternative  ^proach  was  discussed,  namely  design  of  experiment. 


It  was  recently  suggested  to  include  at  least  a  couple  of  benchtop  ejqjeriments  wife  a  very  hi^  concentration 
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Design  of  experiment  (DOE)  is  a  discipline  that  incorporates  statistics  into  the  process  of 
experimentation.  DOE  principles  can  be  used  to  systematically  vary  a  number  of  independent  variables  in 
order  to  evaluate  their  effect  on  dependent  variables  (Hicks,  1982).  In  the  case  of  Orimulsion  behavior  in 
water,  DOE  can  be  used  to  establish  the  minimum  number  of  bench-top  experiments  that  can  be  used  to 
evaluate  the  independent  variables  of  interest  in  a  statistically  valid  manner.  DOE  will  also  enable  one  to 
determine  the  existence  of,  and  even  quantify,  interaction  between  variables. 

3.2  Evaluation  of  Variables 

It  is  important  to  distinguish  between  independent  and  dependent  variables.  Independent  variables  are 
those  parameters  that  will  be  deliberately  controlled  while  dependent  variables  are  those  parameters  that 
are  to  be  measured.  In  preparation  for  this  study  there  were  six  independent  variables  identified  as  having 
a  role  in  bitumen  particle  growth  (see  Section  3.1).  It  is  appropriate  to  ask  at  this  point,  are  these 
independent  variables  appropriate  for  the  purposes  of  this  study?  In  the  Section  3.2.1,  data  and 
observations  from  Battelle’s  previous  experiments  on  Orimulsion  and  those  by  other  researchers  are 
reviewed  with  the  goal  being  able  to  justify  including  all  of  the  variables  selected  for  study. 

3.2.1  Independent  Variables 

3.2.1.1  Water  Salinity  and  Temperature 

Most  oceanic  waters  occur  over  a  relatively  narrow  range  of  salinity  (33-37  ppt).  However,  the  salinity  of 
coastal  waters  can  vary  widely  due  the  competing  influences  of  evaporation,  precipitation  and  freshwater 
river  discharge.  Therefore,  consideration  of  widely  varying  salinity  (0-35  ppt)  is  certainly  justified. 

Temperature  variation  due  to  unequal  heating  with  latitude  and  seasonality  produces  a  significant  range  of 
surface  water  temperatures.  The  ocean  is  warmest  (25-30°C)  in  the  tropical  and  subtropical  regions  and 
coldest  (down  to  -1.7“C)  near  the  poles.  Given  the  latitude  variation  of  the  continental  U.S.,  the  potential 
for  Orimulsion  shipping  into  continental  U.S.  waters  limits  the  temperature  range  to  between 
approximately  10  and  25°C.  Therefore,  our  plan  to  consider  temperature  over  the  range  of  4”C  to  25“C  is 
more  than  adequate. 

The  most  important  influence  of  salinity  and  temperature  on  Orimulsion’s  behavior  is  probably  related  to 
their  effect  on  density.  It  is  well  known  that  the  density  of  seawater  varies  largely  with  salinity  and  is 
further  impacted  by  temperature  (see  Knudson’s  Hydrographic  tables,  e.g.,  in  Riley  and  Chester,  1971). 
Direct  comparison  of  Orimulsion’s  measured  density  at  different  temperatures  (e.g.,  1.0176  at  0°C  to 
1.0048  at  25°C;  Ostazeski  et  al.,  1998)  to  a  seawater  density  table  reveal  that  bitumen  particles  are  more 
dense  than  seawater  with  salinity  <  20  ppt  between  O^C  and  20°C,  and  less  dense  than  seawater  with 
salinity  >  25  ppt  between  25°C  and  30°C.  Thus  considering  no  other  variables,  this  relationship  would 
indicate  that  bitumen  should  sink  when  salinity  <  20  ppt  (between  0-20°C)  and  <  25  ppt  (between  25- 
30“C).  Oppositely,  it  should  float  when  salinity  >  20  ppt  (between  0-20°C)  and  >  25  ppt  (between  25- 
30°C). 

This  float/sink  relationship  was  confirmed  in  the  static  vessel  experiments  by  Febres  et  al.  (1995)  in 
which  bitumen  particles  sank  and  floated  when  salinity  was  below  and  above  20  ppt,  respectively. 
However,  similar  static  vessel  tests  conducted  by  Brown  et  al.  (1995)  at  24.3°C  showed  that  after  an 
initial  period  of  stirring  the  bitumen  floated  in  all  instances  where  the  salinity  exceeded  10  ppt.  Thus, 
there  is  still  some  discrepancy  in  what  to  expect  under  different  salinity  conditions  and  water  salinity  and 
temperature  are  certainly  important  variables  worthy  of  further  study  herein. 
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An  important  effect  of  salinity  may  lie  in  the  influence  that  dissolved  salts  have  on  the  effectiveness  of  the 
surfactant.  This  influence  is  evidenced  by  the  coalescence  of  bitumen  in  seawater  and  lack  of  coalescence 
in  freshwater  (Ostazeski  et  al.,  1997a).  It  can  be  envisioned  that  the  dielectric  potential  of  water  varies 
with  salinity  and  that  charged  ions  in  the  water  (i.e.,  salts)  may  reduce  the  effectiveness  of  the  surfactant 
bound  to  bitumen  particles. 

3.2.1.2  Orimulsion  Concentration 

In  addition  to  the  float/sink  behavior  of  Orimulsion,  there  is  the  issue  of  coalescence  of  bitumen  particles. 
Among  other  factors,  coalescence  must  be  related  to  the  effectiveness  of  the  surfactant  after  being  spilled. 
Orimulsion  originally  contains  approximately  30%  water  (Jokury  et  al.,  1995;  Febres  et  al.,  1995).  When 
it  is  spilt  into  water  the  bitumen  particles  readily  disperse,  in  other  words,  the  water  content  increases 
drastically.  A  substantial  proportion  of  surfactant  is  thought  disassociate  from  the  bitumen  particles  into 
the  water  (Potter,  1995).  This  process  would  seem  to  promote  coalescence  of  bitumen  particles  if  they 
collide.  Theoretically  speaking,  the  collision  frequency  of  bitumen  particles  should  be  a  function  of  (1) 
bitumen  concentration,  (2)  energetics  in  the  system,  and  (3)  the  buoyancy  forces  resulting  from  the 
salinity-temperature  driven  density  differences  between  the  water  and  the  bitumen  (as  described  above). 
Therefore,  our  experiments  must  include  each  of  these  factors.'’ 

Experiments  have  shown  that  when  the  Orimulsion  concentration  in  water  is  high  (>20,0(X)  mg/L)  it 
retains  the  properties  of  the  emulsion  (Febres  et  al.,  1995).  One  can  assume  that  under  such  conditions 
the  surfactant  concentration  in  the  water  and  on  the  bitumen  particle  surfaces  is  high  enough  to  prevent 
coalescence  of  bitumen  particles.  However,  if  the  concentration  is  low  (<10,000  mg/L)  the  Orimulsion 
looses  the  properties  of  an  emulsion  and  bitumen  particles  become  widely  dispersed  (Febres  et  al.,  1995). 
Under  these  conditions  the  surfactant  concentration  in  the  water  is  too  low  to  prevent  coalescence  of  any 
colliding  bitumen  particles  (i.e.,  the  critical  micelle  concentration  is  not  exceeded).  Theoretically 
speaking,  if  the  concentration  of  Orimulsion  spilled  is  very  low  (<1000  mg/L)  the  extremely  low 
concentration  of  bitumen  particles  in  the  water  may  preclude  coalescence  due  to  the  scarcity  of  collisions 
between  particles. 

As  noted  above,  collision  frequency  is  undoubtedly  also  related  to  the  turbulent  energy  of  the  system. 
(This  is  discussed  further  below.)  However,  if  energy  alone  were  responsible  for  collisions  and  the 
resulting  coalescence  of  bitumen  particles,  then  one  would  expect  the  same  degree  of  coalescence  to 
occur  in  freshwater  as  occurs  in  saltwater.  However,  the  flume  experiments  of  Ostazeski  et  al.  (1997a; 
1997b;  1998)  have  shown  that  under  the  same  bitumen  concentration  and  energetic  conditions,  the 
particle  growth  (particularly  surface  patty  formation)  is  much  greater  in  saline  conditions. 

The  reason  for  the  different  responses  in  freshwater  and  saltwater  has  been  theorized  to  be  the  result  of  a 
chemical  ‘denaturing’  of  the  surfactant  in  saltwater  (Ostazeski  et  al.,  1998).  However,  this  theory  remains 
untested  as  there  is  no  evidence  of  the  chemical  change.  It  could  also  be  theorized  that  the  degree  of 
disassociation/partitioning  of  the  surfactant  from  the  bitumen  may  be  greater  in  saltwater  than  in 
freshwater.  However,  Ostazeski  et  al.  (1998)  showed  that  the  dissolved  surfactant  concentration  remained 
the  same  and  stable  throughout  both  the  freshwater  and  saltwater  flume  experiments.  Therefore,  the  third 
factor  controlling  collision  frequency  may  well  lie  in  the  density-float/sink 


However,  based  upon  experiments  involving  clay  particles  in  water,  other  factors  that  undoubtedly  play  at  least 
minor  roles  in  collision  frequency  of  bitumen  particles  are;  (1)  the  particle  distribution  and  shape  (i.e.,  polydispersed 
populations  and  non-spherical  particles  generally  exhibit  greater  collision  potential  than  uniform  particles)  and  (2) 
orthokinetic  coagulation  (i.e.,  the  rate  of  collision  as  smaller  particles  are  ‘captured’  by  rising  or  falling  larger 
particles). 
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relationship  discussed  above.  One  can  envision  how  the  number  of  collisions  might  be  increased  in 
saltwater  where  the  buoyant  bitumen  particles  will  tend  to  rise,  become  more  concentrated  in  the  siuface 
layer,  and  thereby  increase  the  fiequency  of  their  collisions.  On  the  other  hand,  in  fieshwater,  where  the 
bitumen  particles  remain  dispersed  (or  even  sink)  the  colhsion  frequency  is  greatly  reduced.  Perh^s, 
imder  real  spill  conditions  in  freshwater  the  bitumen  particles  would  sink  to  the  bottom  and,  perhrqrs,  even 
coalesced  there  (i.e.,  form  a  sinking  patty).  Not^ly,  in  die  Ostazeski  et  al.  (1997a;  1997b;  1998) 
freshwater  fliune  tank  experiments  (with  a  finite  and  relatively  shallow  water  depth)  the  tank’s  energetics 
kept  the  bitumen  particles  dispersed  and  thereby  may  have  prevented  their  sinking  (and  peih^s  their 
coalescence).  Undoubtedly,  a  range  of  Orimulsion  concentrations  must  be  included  in  our  eiqieriments. 
The  proposed  range  of  17,000  to  1000  mg  Orimulsion/L  water  seems  ^propriate  to  test  the  effect  of 
concentration  on  coalescence.  However,  it  was  recently  suggested  that  a  much  higher  concentration,  ca. 
50,000,  also  be  included  for  observational  purposes  only.  Therefore,  at  least  a  couple  of  tests  wiU  be 
performed  at  this  concentration. 

3.2. 1.3  Particulate  Amount  and  Type 

Particulate  matter  in  oceanic  waters  includes  both  organic  (living  and  dead)  and  inorganic  moieties.  The 
concentration  of  particulates  is  typically  reported  as  the  total  suspended  solids  (TSS).  For  our  purposes, 
the  importance  of  TSS  in  the  water  into  which  Orimulsion  is  spilled  probably  lies  with  their  potential  to 
increase  coalescence  of  bitumen  particles  by  increasing  the  potential  for  collisions.  This  phenomenon 
was  observed  by  Brown  et  al.  (1995)  who  showed  fliat  when  excessively  large  concentrations  of  (Tampa 
Bay)  sediment  was  added  to  agitated  vessels  significant  quantities  of  bitumen  became  adhered  to  the  filner 
particles  and  rapidly  sank  (due  to  the  relatively  high  specific  gravity  of  tbe  mineral-bitumen  particles).  It 
was  envisioned  that  suspended  sediments  in  water  provided  a  surfactant-fiee  surface  to  which  bitumen 
particles  could  adhere  upon  a  collision. 

The  concentration  of  TSS  in  water  will  be  a  function  of  the  particle  size  and  the  turbulent  energy  in  the 
system  (i.e.,  energetics).  In  addition,  in  the  case  of  fieshwater  the  concentration  of  some  particulates  will 
also  be  dependent  upon  their  solubility  (e.g.,  calcite  added  to  fieshwater  will  be  subject  to  some 
dissolution).  Furthermore,  since  energetic  conditions  in  our  experiments  will  vary  we  can  expect  the  TSS 
concentration  to  vary.  Therefore,  we  propose  to  add  a  given  amormt  of  particulate  matter  to  each 
experiment  that  will  represent  the  ‘high’  particulate  load.  The  Tow’  load  amount  to  be  added  will  be 
zero.  The  combined  effect  of  the  energetics  and  particle  size  of  this  material  will  then  determine  the 
amount  of  particulate  in  suspension.  What  amormt  is  appropriate  to  add? 

Open  ocean  concentrations  of  TSS  are  generally  in  the  0.5-250  pg/L  range  (Riley  and  Chester,  1971). 
Concentrations  in  areas  of  deltaic  sedimentation  are  much  higher  (e.g.,  Mississippi  River  average  270 
mg/L;  Milliman,  1990).  We  propose  to  select  a  ‘high’  concentration  of  arormd  270  mg/L.'"  We  envision 
only  a  certain  amormt  of  this  will  stay  in  suspension,  especially  rmder  the  low  energetic  conditions. 

What  kind  of  material  should  be  added?  Theoretically  speaking,  smaller-sized  particulate  matter  would 
seem  to  favor  collisions  and  adherence  of  bitumen  due  to  the  increased  sru&ce  area  and  ease  of 
suspension  of  smaller  particulates,  e.g.,  clay-  and  silt-sized  particles  (<  2  and  2-62.5  pm,  respectively). 
This  was  also  confirmed  in  the  study  by  Brown  et  al.  (1995)  who  showed  that  these  fine  grained 
particulates  (<  75  pm)  demonstrated  ^proximately  three  orders  of  magnitude  greater  adhesion  potential 
than  coarser  grained  particulates  (>  75  pm).  Small  particulates  also  comprise  the  majority  of  the 
lithogenous  sediment  load  carried  to  coastal  waters  by  rivers  and  these  particles  also  have  the  longest 
residence  time  in  seawater.  Thus,  only  clay-to-silt  size  particles  will  be  considered  in  these  experiments. 

'  Note  that  the  35,710  to  142,860  mg/L  sediment  concentrations  used  in  the  experiments  by  Brown  et  al.  (1995)  far 
exceeded  the  typical  suspended  solid  (TSS)  load  in  Tampa  Bay  waters  which  ranges  from  20-120  mg/L. 
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What  kind  of  clay-silt  size  particles  will  be  used?  It  had  been  suggested  that  two  types  of  particulate 
material  be  used  in  these  ejq)eriments,  kaolinite  and  calcite.  Kaolinite’s  selection  is  reasonable,  as  it 
comprises  a  significant  proportion  (10-20%)  of  the  clay-sized  materials  in  U.S.  coastal  marine  sediments. 
The  impact  of  additional  particles  (other  than  bitumen)  in  our  experiments  will  confound  the 
measurement  of  bitumen  particle  size.  Kaolinite,  like  all  clays,  are  platy  minerals  whose  crystal  structure 
consists  of  layers  or  sheets.  Thus,  the  size  of  an  individual  kaolinite  particle  varies  significantly 
depending  on  the  orientation  of  the  plates  and  the  number  of  plates  in  a  ‘stack’.  Our  preliminary 
investigations  have  revealed  that  the  particle  size  distribution  of  a  kaolinite  slurry  exhibits  maxima  around 
2  |im  and  25  pm.  We  believe  these  maxima  represent  the  ‘side’  and  ‘top’  views  of  kaolinite  particles. 
These  populations  of  particles  caimot  be  distinguished  firom  bitumen;  therefore,  it  must  be  realized  that 
for  the  experiments  in  which  kaolinite  is  added  the  particle  diameters  being  measured  will  be  a  composite 
of  bitumen,  kaolinite,  and  potentially  bitumen/kaolinite  particles. 

It  had  also  been  suggested  that  we  consider  the  influence  of  a  non-lithogenous  particulate  material,  e.g.,  a 
precipitated  material.  Calcite  is  fire  most  common  carbonate  mineral  which  forms  fix>m  the  precipitation 
of  CaCOa  fix)m  seawater.  Since  calcite  is  not  carried  by  rivers  to  larger  water  bodies  and  deposited,  it 
t3^ically  does  not  occur  as  TSS  in  most  settings.  In  areas  where  carbonate  sedimentation  predominates 
(e.g.,  coastal  Florida)  it  is  possible  to  get  calcite  particles  into  suspension  through  disturbing  of  the  bottom 
sediments.  Therefore,  we  will  still  consider  including  calcite.  However,  one  important  foctor  that  has  not 
been  considered  in  our  earlier  discussions  is  the  fact  that  if  we  were  to  use  a  second  type  of  particulate 
material,  the  particle  size  of  it  would  be  different  firom  the  other  (e.g.,  kaolinite).  Calcite  will  also 
partially  dissolve  in  our  fineshwater  experiments.  Since  particle  size  measurements  are  already  going  to  be 
confounded  by  the  presence  of  kaolinite,  we  will  further  confound  the  particle  size  results  by  using  a 
second  type  of  particulate  material.  Nonetheless,  it  is  considered  valuable  to  include  calcite  in  order  to 
address  the  issue  of  whether  or  not  the  type  of  particulate  matter  may  influence  the  coalescence 
phenomenon.  As  will  be  evident  in  the  Section  3.3,  the  use  of  both  kaolinite  and  calcite  introduces  a 
discrete  variable  into  the  experimental  design.  In  doing  so,  the  options  for  the  design  become  limited, 
therefore,  we  expect  that  the  addition  of  calcite  will  be  limited  to  the  screening  phase  of  the  study. 

3.2. 1.4  Energetics 

Previous  benchtop  experiments  involving  Orimulsion  have  been  performed  under  static  conditions  or 
under  gentle  agitation  by  magnetic  stirrers.  No  attempt  was  made  to  quantify  the  energy  in  the  latter 
systems.  In  our  study  we  intend  to  introduce  some  degree  of  turbulence  to  all  of  our  benchtop 
experiments  (i.e.,  no  static  tests).  The  objective  of  evaluating  this  variable  is  to  determine  the  impact  that 
turbulent  energy  has  on  the  frequency  of  collisions  between  bitumen  particles  and  between  bitumen  and 
clay  particles.  Therefore,  the  method  used  need  only  increase  the  propensity  for  collisions. 

Naturally,  there  are  a  variety  of  ways  in  which  turbulence  can  be  introduced.  We  have  considered  many 
different  options.  We  intentionally  have  not  considered  any  method  that  reqirired  putting  some  t5q)e  of 
agitator  into  the  vessels  (e.g.,  oscillating  rings,  propellers,  etc.)  due  to  the  ‘stickiness’  of  Orimulsion.  We 
have  also  considered  rotating  flasks  and  based  on  earlier  experience  have  eliminated  this  option. 
(Previous  experiments  with  Orimulsion  in  seawater  have  produced  a  bitumen  ‘plug’  in  these  vessels.)  We 
have  also  had  to  consider  the  need  to  conduct  a  fairly  large  number  of  experiments  on  moderate  volumes 
of  water.  Therefore,  fiom  a  practical  standpoint  we  propose  to  introduce  turbulence  to  our  benchtop 
experiments  by  rtsing  an  orbital  shaker  table  apparatus  (described  in  Section  4.1).  A  brief  discussion  here 
is  directed  at  understanding  the  turbulence  within  such  a  system. 

A  turbulence  spectrum  will  include  eddies  of  different  amounts  of  energy  and  sizes.  According  to 
turbulence  theory  the  largest  eddies  will  cascade  into  smaller  eddies  and  the  smallest  will  dissipate  by 
being  converted  into  heat  by  the  viscosity  of  the  water  (Delvigne,  1993).  In  our  proposed  shaker  table 
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experiments  the  laiger  eddies  created  will  be  of  the  same  order  of  magnitude  as  the  disturbance,  i.e.,  the 
diameter  of  the  container.  Recently,  Intervep  researchers  (Masciangioli  et  al.,  personal  communication, 
1998)  have  conducted  orbital  shaker  table  studies  in  which  the  cascade  process  is  believed  to  occur  and 
the  kinetic  energy  (K)  of  the  water  due  to  its  rotation  (i.e.,  not  its  oscillation  due  to  the  surfece  wave) 
could  be  estimated  by  the  following  equation: 

K  =  7tpa)^HRV4 

Where  p=density  of  fluid,  a)=mean  angular  velocity  of  the  fluid,  H=water  depth  and  R=radius  of  vessel. 
They  demonstrated  that  in  a  IL  vessel  with  H=12.7  and  R=5  cm,  rotated  at  120  tpm,  die  water  had  an 
oscillation  firequency  of  25  rpm  (reduced  fiom  the  table’s  rpm’s  due  to  inertia,  friction  and  other  effects) 
and  the  K  could  be  calculated  fiom  the  equation  above  to  be  0.52  Joules/m^.  The  mean  tangential 
velocity  of  a  water  particle  5  cm  fiom  the  center  of  the  vessel  was  determined  to  be  13  cm/sec.  This  was 
considered  comparable  to  a  1.44  cm  amplitude  progressive  wave  having  an  orbital  velocity  of  16.5 
cm/sec. 

We  have  conducted  some  preliminary  experiments  using  our  prototype  benchtop  qiparatus  (described  in 
Section  4.1).  We  propose  that  the  following  operating  conditions  be  used:  diameter  of  oscillation  of  1.9 
cm,  shaker  table  fiequencies  of  10  and  100  rpm,  shaker  table  oscillation  periods  of  6.0  and  0.60  sec, 
radius  of  3L  water  of  9  cm,  height  of  3L  of  water  is  14  cm.  Calculations  like  those  described  above  were 
performed  on  these  operating  conditions.  These  resulted  in  a  K  of  1.096  J/m^.  The  mean  tangential 
velocity  of  a  water  particle  8  cm  fiom  the  center  of  the  vessel  was  determined  to  be  17.5  cm/sec.  This 
was  considered  comparable  to  a  2.1  cm  amplitude  progressive  wave  having  an  orbital  velocity  of  12.0 
cm/sec. 

3.2. 1. 5  Summary  of  Independent  Variables  to  be  used  in  Benchtop  Experiments 

Table  3-1  summarizes  the  evaluation  of  the  independent  variables  that  we  propose  to  include  in  our 
benchtop  experiments.  All  but  one  of  the  variables  to  be  included  are  continuous  variables,  i.e.,  they  can 
have  any  value  over  some  numerical  range.  The  complicating  factors  introduced  by  the  inclusion  of  the 
discrete  variable,  particulate  type,  have  been  discussed.  It  was  nonetheless  decided  to  include  this 
variable  in  the  screening  phase  of  our  study  (See  below). 


Table  3-1:  Summary  of  Independent  Variables  to  be  use  in  our  Benchtop  Experiments 


EXPERIMENTAL  RANGE 

INDEPENDENT  VARIABLES 

FACTOR 

TYPE 

LOW 

fflGH 

CENTER 

UNITS 

Salinity 

continuous 

0 

35 

ppt 

Water  temperature 

continuous 

4 

25 

14.5 

"C 

Orimulsion  concentration 

continuous 

1000 

17000 

Particulate  loading 

continuous 

0 

270 

135 

mg/L 

Particulate  Type 

discrete 

kaolinite 

calcite 

— 

— 

Water  energetics 

continuous 

10 

100 

45 

rpm 

3.2.2  Dependent  Variables 
3. 2. 2. 1  Bitumen  Particle  Size 

The  primary  dependent  variable  that  was  to  be  measured  is  bitumen  particle  size.  The  method  for 
collection  and  measurement  of  bitumen  from  the  benchtop  vessels  is  described  in  Section  5.2.  As 
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mentioned  above,  bitumen  particle  size  measurements  will  be  complicated  by  the  addition  of  kaolinite 
and  calcite  in  some  of  the  experiments.  Our  preliminary  experiments  have  shown  that  it  will  be  extremely 
difficult  to  describe  bitumen  particles  separately  from  kaolinite  or  calcite,  therefore,  any  particle  size 
measurements  will  be  a  composite  of  bitumen,  kaolinite  (or  calcite),  and  bitumenykaolinite  (or 
bitumen/calcite)  particles.  Nonetheless,  with  careful  interpretation  and  knowledge  of  kaolinhe’s  and 
caldte’s  particle  size,  it  should  still  be  possible  to  qualitatively  recognize  changes  in  overall  particle  size. 
However,  absolute  bitumen  particle  size  carmot  be  determined  in  the  experiments  that  will  include 
kaolinite.  Those  experiments  that  will  not  include  any  kaolinite  or  calcite  will  allow  for  bitumen  particle 
growth  to  be  quantified. 

The  need  to  measure  changes  in  particle  size  over  the  course  of  the  experiments  requires  that  either  the 
measurements  be  performed  at  discrete  intervals  throughout  the  experiments  or  that  experiment  duration 
become  another  independent  variable.  The  former  seems  tike  a  more  cost-effective  ^proach  since  far 
fewer  experiments  will  need  to  be  conducted.  We  propose  to  measure  particle  size  (bitumen,  clay/calcite, 
or  bitumen/clay/calcite)  of  the  suspended  material  in  intervals  of  1,  3,  8  and  24  hours  following  the 
release  of  Orimulsion  in  our  benchtop  studies.  These  four  intervals  should  provide  the  information 
necessary  to  determine  if  particle  growth  is  occurring  and  at  what  relative  rate  for  the  variable  conditions 
of  our  experiments.  One  potential  problem  with  this  tqrproach  is  that  after  24  hours  it  is  possible  that 
there  will  be  too  few  suspended  particles  to  provide  a  reliable  particle  size  measurement. 

Another  confounding  aspect  of  particle  size  measurements  will  be  the  formation  of  a  wide  range  of 
particle  sizes  (e.g.,  fiom  2  micron  clay  particles  to  10  cm  diameter  patties).  Measurement  of  particles 
over  such  a  wide  dynamic  range  is  virtually  imposable.  Only  via  techniques  such  as  microscopy  or 
sedimentation  can  this  wide  range  of  particles  be  measured.  However,  these  techniques  carmot  be  used  in 
the  case  of  bitumen  since  they  would  undoubtedly  change  the  particles’  environment  to  the  point  where 
the  particle  size  would  continue  to  change.  Thus,  the  measurement  technique  must  be  relatively  fest.  The 
most  qjpropriate  techruque  is  laser  diffiaction  but  this  technique  is  unable  to  ‘see’  particles  over  such  a 
wide  dynamic  range.  (Even  the  most  expensive  imits  can  only  see  fiom  0.4  to  2000  microns).  The  widest 
dynamic  range  available  on  oiu  instrument  is  1.5  and  600  microns  (see  Section  5.1.4).  Our  experience 
tells  us  that  particles  well  above  600  microns  can  be  formed  in  these  experiments.  Nonetheless,  we  will 
measure  particle  size  throughout  our  experiments  with  the  intention  being  to  recognize  changes  rather 
than  absolute  particle  diameters. 

3.2.2.2  Mass  Balance  Measurements 

We  believe  that  bitumen  particle  size  should  not  be  the  only  dependent  variable  to  be  measured  in  these 
experiments.  A  variety  of  chemical  and  physical  measurements  can  be  made  that  will  provide  additional 
insight  into  the  fate  of  Orimulsion  when  spilled  under  a  variety  of  conditions.  After  all,  this  is  the 
primary  goal  of  these  studies.  First  and  foremost  among  these  is  an  attempt  to  ‘mass  balance’  the  end 
products. 

A  mass  balance  of  bitumen’s  fate  was  attempted  in  the  previous  flume  experiments  conducted  in 
Battelle’s  4000L  tank.  In  the  seawater  study  it  was  estimated  that  at  the  end  of  the  168  hour  test  that  65% 
of  the  bitumen  occurred  in  the  surface  patty,  ~35%  was  adhered  to  the  walls  near  the  water  level,  and 
~1%  was  still  within  the  water  colurtm.  We  intend  to  attempt  a  more  thorough  reconciliation  of  the  fete 
of  Orimulsion  in  our  benchtop  experiments.  Our  preliminary  experiments  using  the  prototype  apparatus 
have  shown  that  Orimulsion’s  bitumen  will  either 

(1)  form  a  floating  patty, 

(2)  form  a  sinking  patty. 
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(3)  coalesce  with  kaolinite  and  sink  to  fonn  non- 
cohesive  grains  of  bitumen/kaolinite, 

(4)  dissolve  into  the  water, 

(5)  remain  suspended  in  the  water,  or 

(6)  stick  to  the  vessel  walls  in  the  vicinity  of  the 
water  line  (i.e.,  form  a  “bathtub  ring”). 

We  propose  to  quantify  the  bitumen  that  ends  up  in  each  of  these  six  ‘compartments’  (see  Section  5.2.6). 

The  weight  of  the  floating  patty  will  be  measured  at  the  conclusion  of  the  experiment.  It’s  moisture 
content  will  be  determined  in  order  to  ‘back  out’  the  mass  of  bitumen  in  the  patty.  If  a  sinking  patty  has 
formed  it  will  be  weighed  and  analyzed  in  the  same  maimer.  Any  sediment  that  accumulates  in  the 
bottom  of  the  vessels  after  24  hours  will  be  collected  and  analyzed  to  determine  the  concentration  of  total 
petroleum  hydrocarbons  (TPH)  that  occurs  in  the  ‘sink’  fiaction.  The  TPH  fliat  occmrs  in  suspension  will 
be  determined  after  1  and  8  hours.  The  concentration  of  TPH  that  is  dissolved  (filtered  water)  and 
remains  suspended  (filtered  material)  in  the  water  column  after  24  hours  will  be  determined  on  an  aliquot 
of  water.  Finally,  the  amount  of  bitumen  that  has  adhered  to  flie  vessel  will  be  determined 
gravimetrically. 

The  feet  that  some  bitumen  will  be  removed  fijom  the  vessels  in  the  course  of  the  experiments  for  particle 
size  and  TPH  measurements  implies  that  some  small  percentage  of  the  bitumen  will  not  be  accounted  for 
in  the  mass  balance.  Nonetheless,  the  results  will  provide  a  semi-quantitative  and  relative  sense  of  just 
where  the  oiganic  matter  is  going.  These  semi-quantitative  results  will  provide  the  basic  information 
necessary  to  predict  Orimulsion’s  behavior  under  the  wide  range  of  conditions  to  be  studied  and  also  help 
to  focus  the  meso-scale  experiments.  We  believe  that  herein  may  lie  the  greatest  practical  yield  of 
information  from  the  benchtop  experiments  proposed  herein. 

3.2.23  Fate  of  Polycyclic  Aromatic  Hydrocarbons 

Polycyclic  aromatic  hydrocarbons  are  perhaps  the  most  important  chemicals  of  concern  in  terms  of  their 
potential  impact  to  marine  and  coastal  biota.  These  experiments  provide  the  opportunity  to  investigate  the 
influence  of  several  variables  on  determining  where  the  PAH’s  would  end  up  in  the  event  of  an 
Orimulsion  spill.  We  propose  to  determine  the  concentration  of  PAH  that  occur  dissolved  in  the  water 
and  those  which  occur  in  the  sediment  at  the  conclusion  of  the  benchtop  experiments. 

3.2.2.4  Summary  of  Dependent  Variables  to  be  Measured  in  Benchtop  Experiments 

In  summary  of  the  discussion  above.  Table  3-2  lists  the  dependent  variables  that  will  be  measured  in  the 
course  of  the  benchtop  experiments. 
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Table  3-2:  Inventory  of  proposed  independent  and  d^endent  variables. 


DEarail^ENT  VAittABI^ 

particle  size  in  suspension  after  1  hour  -  pm 

particle  size  in  suspension  after  3  hours  -  pm 

particle  size  in  suspension  after  8  hours  -  pm 

particle  size  in  suspension  after  24  hours  -  pm 

Mass  of  floating  bitumen/patty  after  24  hours  -  g 

Mass  of  sinking  bitumen/patty  after  24  hours  -  g 

Mass  of  bitumen  adhered  to  glass  walls  after  24  hours  -  g 

TPH  in  sediment  after  24  hours  -  mg/kg 

TPH  in  suspension  after  1  hours  -  mg/L 

TPH  in  suspension  after  8  hours  -  mg/L 

TPH  in  suspension  after  24  hours  -  mg/L 

Dissolved  TPH  in  water  after  24  hours  -  pg/L 

Dissolved  PAH  in  water  after  24  hours  -  pg/L 

3.3  Design  of  Experiment  Principles 

There  are  three  basic  levels  to  DOE,  (1)  screening  level,  (2)  two-level  &ctorial  design  and  (3)  three-level 
response  surface  design.  The  screening  level  experiments  are  intended  to  identify  those  variables  that 
have  the  greatest  ejffect  on  the  outcome  and  hereby  reduce  the  number  of  variables  for  a  more 
sophisticated  two-  or  three-level  design  study.  In  our  case,  screening  would  help  identify  which  of  the 
five  independent  variables  (Table  3-1)  have  die  greatest  effect  on  the  dependent  variables  (Table  3-2). 
Because  we  are  monitoring  multiple  dependent  variables  some  judgement  may  have  to  be  made  regarding 
the  ‘most  important’  ones  if  it  is  foimd  that  they  are  controlled  by  different  independent  variables.  In  this 
study  it  will  be  pradent  to  perform  a  set  screening  level  experiments  in  order  to  potentially  reduce  the 
number  of  independent  variables  in  some  additional  two-  or  three-level  experimental  designs. 

3.3.1  Screening  Designs 

There  are  two  types  of  screening  techniques  Plackett-Burman  (PB)  and  fiactional  factorial  (FF),  with  the 
former  being  much  more  widely  used  do  to  its  ease  of  interpretation.  PB  designs  are  multi-variable,  two- 
level  designs  that  can  be  used  to  evaluate  any  number  of  variables  up  to  four  less  than  the  number  of 
experiments.  Given  the  six  independent  variables  determined  for  these  experiments  (Table  3-1)  we 
propose  to  screen  using  a  12-run  PB  design  as  shown  in  Table  3-3.  In  Table  3-3  the  ‘+’  and  signs 
represent  the  high  and  low  values  for  each  of  the  dependent  variables.  Assigning  these  values  to  our  list 
of  five  variables  yields  the  design  shown  in  Table  3-4. 
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Table  3-3;  Twelve-run  Plackett-Burman  design 


XIO  Xll 


Table  3-4:  Factors  assigned  to  the  first  six  columns  of  Table  3-3  and  sorted  by  temperature. 


In  Table  3-4  the  first  six  columns  contain  the  assigned  ‘extreme’  values  for  the  independent  variables. 
(The  extreme  Orimulsion  concentration,  50,000  mg/L,  is  not  included  in  this  design-see  below).  Using 
these  extremes  is  considered  to  be  an  important  part  of  any  DOE  as  it  will  make  it  easier  to  recognize  the 
variables  which  have  the  greatest  effect.  Examination  of  Table  3-3  reveals  that  there  is  a  balance  among 
all  of  the  variables  so  that  each  variable  is  paired  with  each  other  variable’s  extreme  values.  While  it 
would  be  most  appropriate  to  conduct  the  experiments  (trials)  sequentially  and  in  a  random  order  this  is 
considered  impractical  given  the  need  to  minimize  costs.  Since  temperature  is  the  most  difficult 
parameter  to  control,  we  propose  to  perform  two  sets  of  6  experiments  to  be  conducted  at  4°C  and  25°C. 

The  last  five  columns  in  Table  3-4,  the  unassigned  columns,  will  be  used  to  measure  the  “noise”  (or 
standard  error  of  the  fector  effect;  Sfe)  within  the  system.  The  effect  of  these  unassigned  factors  will  also 
provide  a  means  of  assessing  the  presence  of  a  large  interaction  between  variables. 

The  dependent  variables  that  will  be  quantified  in  the  course  of  the  12  PB  design  experiments  include  all 
of  those  listed  in  Table  3-2.  Calculation  of  the  t-value  and  95%  confidence  limits  of  Ae  Sfe  will  provide  a 
means  of  identifying  and  ranking  tiiose  independent  variables  that  have  a  significant  effect  on  each  of  the 
dependent  variables  being  measured.  Interpreting  which  independent  variables  lead  to  the  highest  or 
lowest  concentrations  of  TPH  or  PAH,  or  greatest  patty  weight  after  24  hours  will  be  straightforward. 
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However,  an  evaluation  of  the  particle  growth,  i.e.,  particle  size  at  different  times  in  the  course  of  the 
experiments  (1,  3,  8  and  24  hours)  wiU  be  more  difficult  since  it  requires  qualitative  analysis  of  four 
different  results.  We  intend  to  synthesize  the  particle  size  results  for  each  individual  time  period  in  order 
to  identify/rank  the  most  important  variables  leading  to  increased  particle  size  throughout  the  entire 
experiment.  Interpretation  of  diese  variables  on  particle  growth  will  be  further  complicated  by  the 
presence  ofkaolinite/calcite  in  some  of  these  runs. 

It  can  be  anticipated  that  the  PB  screening  process  will  help  to  identify  the  three  to  five  most  significant 
independent  variables  that  effect  each  of  the  dependent  variables  over  24  hours.  However,  as  mentioned 
above,  because  we  are  monitoring  multiple  dependent  variables  some  judgement  will  have  to  be  made 
regarding  the  ‘most  important’  dependent  variables,  particularly  if  it  is  found  that  they  are  controlled  by 
different  independent  variables.  In  any  case,  in  having  eliminated  some  variables  fiom  further 
consideration,  we  have  two  options  by  which  to  proceed,  namely  a  two-level  fiictorial  design  experiment 
or  a  three-level  response  surfiice  design  (i.e.,  a  face-centered  cube  design  or  a  Box-Behnken  spherical 
design).  These  two  scenarios  are  described  below  for  completeness.  However,  for  reasons  described  at 
the  beginning  of  Section  3.3.3,  we  believe  that  we  will  perform  only  tire  three-level  response  surfece 
design. 

In  addition  to  the  PB  screening  design  shown  in  Table  3-4,  it  is  proposed  that  four  additional  experiments 
be  conducted  at  extremely  high  concentrations  of  Orimulsion,  namely  50,000  mg/L.  Four  experiments 
are  proposed  in  which  the  primary  independent  variables  to  be  investigated  are  salinity,  temperature  and 
energetics  (Table  3-4’).  The  dependent  variable  results  of  these  investigation  will  not  be  include  in  the 
statistical  results  of  the  screening  study  hut  rather  they  will  be  investigated  and  described  qualitatively. 
Experiments  at  this  high  concentration  of  Orimulsion  are  intended  to  provide  some  basis  upon  which  to 
predict  its  response  in  the  immediate  vicinity  of  a  release  or  spill,  i.e.,  where  concentrations  would  be 
their  highest. 

Table  3-4’:  High  Orimulsion  concentration  experiments  to  be  conducted  concurrently  with  those  in 

Table  3-4. 


Trial  Dependent  S  Temp.  Orim.  cone.  Part  Loadmg  Energetics  Part  Type 
Variable  (ppt)  (”C)  (mg/L)  (mg/L)  (tpm) 

13 

_ 1 

35 

4 

50,000 

0 

10 

— 

14 

0 

4 

50,000 

0 

100 

~ 

15 

35 

25 

50,000 

0 

100 

— 

16 

0 

25 

50,000 

0 

10 

— 

3.3.2  Two-level  Factorial  Design 

A  two-level  factorial  design  approach  is  used  to  estimate  all  factor  effects  in  linear  interaction  effects. 
(This  qjproach  must  be  used  if  we  were  to  change  our  approach  to  include  any  discrete  variables,  e.g., 
particulate  type.)  The  variables  with  significant  effects  identified  by  the  PB  screening  experiments  will 
be  used  to  constract  either  a  three-,  four-,  five-  or  six-fkrtor  two-level  factorial  design  leading  to  a  total  of 
8,  16,  32  or  64  additional  experiments.  We  fully  anticipate  that  the  PB  screening  can  be  used  to  limit  the 
number  of  significant  variables  to  three  or  four,  thus  limiting  the  number  of  additional  experiments  to  8  or 
16.  As  an  example  to  this  approach.  Table  3-5  shows  the  four-factor  two-level  factorial  design  assuming 
salinity,  temperature,  Orimulsion  concentration,  and  energetics  were  identified  as  being  significant  in  the 
PB  screening. 
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Table  3-5:  Example  of  a  four-factor  two  level  factorial  design. 


XI  X2  X3  XI  X2  X3  X4 

Run# 

m 

m 

Salinity 

Gone. 

Energetic 

s 

Batch 

1 

-1 

-1 

-1 

-1 

0 

4 

1000 

10 

1 

2 

1 

-1 

-1 

-1 

35 

4 

1 

3 

-1 

1 

-1 

-1 

0 

25 

1000 

10 

2 

4 

1 

1 

-1 

-1 

35 

25 

1000 

10 

2 

5 

-1 

-1 

1 

-1 

0 

4 

17000 

10 

1 

6 

1 

-1 

1 

-1 

35 

4 

17000 

10 

1 

7 

-1 

1 

1 

-1 

17000 

2 

8 

1 

1 

1 

-1 

35 

25 

2 

9 

-1 

-1 

-1 

1 

0 

4 

1000 

100 

1 

10 

1 

-1 

-1 

1 

35 

4 

1000 

100 

1 

11 

-1 

1 

-1 

1 

0 

25 

1000 

100 

2 

mm 

1 

1 

-1 

1 

35 

25 

2 

13 

-1 

-1 

1 

1 

0 

4 

17000 

100 

1 

14 

1 

-1 

1 

1 

35 

4 

17000 

1 

-1 

1 

1 

1 

25 

17000 

100 

2 

16 

1 

1 

1 

1 

35 

25 

17000 

100 

2 

No  blocking  or  randomization  would  be  used  but  rather  the  experiments  would  be  conducted  in  two 
‘batches’  depending  on  the  temperature  (4  or  25°C),  again,  in  order  to  reduce  costs.  As  was  the  case  in 
the  PB  screening  the  same  set  of  dependent  variables  (Table  3-2)  would  be  measured  for  each  experiment. 
The  fector  effect  and  interaction  effects  of  these  four  (or  other  number  of)  factors  will  be  ^culated 
according  to  standard  methods  that  would  produce  a  measure  of  the  effect.  The  unassigned  factor  effects 
would  be  used  to  calculate  random  error  in  the  system  and  95%  confidence  limits.  Evaluation  of  the 
results  would  proceed  as  described  above  for  the  PB  screening.  If  the  two-level  fectorial  design  is 
performed,  the  most  important  independent  variables  identified  would  be  used  to  select  a  set  of  conditions 
that  would  be  used  in  the  meso-scale  flume  experiments  (see  below). 


3.3.3  Three-level  Response  Surface  Design 
3.3.3. 1  Overview  and  Justification  for  this  Approach 

A  three-level  response  surface  design  is  preferred  over  the  two-level  factorial  design  because  the  former 
will  allow  for  the  recognition  and  quantification  of  the  non-linear  (i.e.,  curvature)  effects.  In  other  words, 
two-level  factorial  design  can  only  recognize  changes  over  a  linear  response  curve.  However,  our 
preliminary  experiments  and  previous  work  by  others  give  us  reason  to  believe  that  the  properties  of 
Orimulsion  will  cause  it  to  not  necessarily  respond  linearly.  For  example,  at  some  salinity  value  between 
32  and  4  7oo  it  might  be  anticipated  that  particle  growth  and  patty  formation  will  occur  but  at  only  a 
slightly  lower  salinity  particle  growth  and  patty  formation  will  not  occur.  Thus,  the  response  may  be  non¬ 
linear. 

In  order  recognize  this  and  other  non-linear  effects  there  can  be  no  discrete  variables  among  the 
independent  variables  being  evaluated  in  a  three-level  response  surface  design.  (This,  therefore,  is 
another  good  reason  to  not  complicate  the  study  by  including  both  kaolinite  and  calcite.)  This  is  because 
if  curvature  effects  are  to  be  estimated,  each  independent  variable  must  have  at  least  three  levels.  In  our 
opinion,  there  is  no  realistic  intermediate  level  in  Ae  case  of  particulate  type  (i.e.,  50:50  kaolinite:  calcite 


Appendix  I 
Page:  14 


is  not  a  real  world  situation).  It  is  fiar  these  reasons  that  we  propose  to  perform  the  three-level  response 
surface  design  on  the  most  critical  continuous  independent  variables  identified  in  the  PB  screening. 

As  a  result  of  considering  intermediate  values  for  continuous  variables,  three-level  response  sur&ce 
designs  are  better  able  to  describe  the  responses  throughout  the  entire  eTqrerimental  ranges.  Hus  is 
accomplished  by  the  use  of  quadratic  polynomials  to  match  the  response  surfaces.  Ihis  will  permit  more 
accurate  prediction  of  the  responses  even  at  factor  combinations  not  actually  nm.  For  example  i^  as  we 
believe,  salinity  will  play  an  important  role  in  particle  growth  or  floating  patty  weight,  the  response 
siufece  design  will  permit  us  to  identify  the  ‘key’  salinity  value  at  which  these  responses  are  most 
effected.  Furthermore,  this  approach  will  allow  us  to  narrow  the  range  over  which  each  variable  has  a 
sigiuficant  effect  on  the  dependent  variables  of  interest.  These  results  will  permit  an  identification  and 
narrowing  of  the  conditions  to  be  used  during  the  subsequent  meso-scale  studies  of  Phase  III. 


3.3.3.2  Types  of  Three-level  Response  Surface  Designs 

There  are  two  types  of  three-level  response  siufece  designs,  namely  a  fece-centered  cube  design  or  a  Box- 
Behnken  spherical  design.  The  former  considers  a  cubical  experimental  region  and  the  latter  considers  a 
spherical  experimental  region.  The  best  way  to  consider  this  is  to  imagine  a  three  variable  problem.  The 
shape  describing  the  range  for  each  variable  would  be  a  cube,  with  the  high  and  low  values  of  each  of  the 
three  variables  occurring  at  the  comers  of  the  cube.  This  cube  would  have  8  comers,  6  &ce  centers,  and 
one  center  point.  This  design  would  include  all  values  within  the  ranges  for  each  variable  and  is  tire  basis 
for  the  fece-center  cube  design.  Connecting  only  the  midpoints  of  these  same  three  variable  ranges  would 
form  a  polyhedron  that  r^proximates  a  sphere  (within  the  original  cube).  This  polyhedron  would  have  12 
edge  points  and  one  center  point.  This  design  would  exclude  the  comers  of  the  cube,  i.e.,  the  high  and 
low  values  for  each  variable  and  is  the  basis  for  the  Box-Behnken  spherical  design. 

The  fece-centered  cube  design  is  considered  more  ‘bold’  than  the  Box-Behnken  design.  The  fece- 
centered  cube  design  may  be  more  qrpropriate  for  the  problem  at  hand  since  it  includes  tire  entire  range 
for  each  independent  variable.  It  seems  pmdent  to  not  ignore  the  ‘extreme’  conditions  given  that  they 
may  occur  in  the  ‘real  world’.  For  example,  it  is  possible  that  a  situation  exists  were  Orimulsion  is 
released  into  35  7oo  seawater,  at  25'’C,  with  no  inorganic  particulates,  and  ‘mild’  energetics.  Therefore, 
we  recommend  that  a  fiice-centered  cube  design  be  used  in  this  part  of  Phase  H.  The  number  of 
sigrtificant  variables  identified  in  the  PB  screening  will  determine  the  number  of  experiments  necessary  in 
a  face-centered  cube  design  according  to  Table  3-6. 

Table  3-6:  Size  of  response  surface  model  necessary  for  2  to  5  independent  variables. 


Number  of  Experiments  required - 

Variables  Face-centered  cube 

design 

2 

11 

3 

17 

4 

31 

5 

45 

If  all  five  of  the  independent  variables  (Table  3-1)  are  identified  as  important  then  45  additional 
experiments  will  be  necessary.  If  any  can  be  eliminated  then  31  or  fewer  experiments  will  be  necessary. 
This  number  will  still  provide  a  modest  niunber  of  degrees  of  freedom  for  calculating  error. 
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Let  US  assume  that  the  PB  screening  has  determined  that  only  salinity,  temperature,  and  bitumen 
concentration  are  important  in  determining  most  of  the  dependent  variables  (Table  3-2).  Evaluating  these 
three  independent  variables  would  require  17  experiments  performed  according  to  the  three-factor  face- 
centered  cube  design  shown  in  Table  3-7. 


Table  3-7:  Hypothetical  three-factor  face-centered  cube  design. 


Btih# 

fSIISIIill 

ISI 

Salinity 

IBS 

Gone. 

Batch 

agam 

mg/L 

1 

Comer 

- 

- 

- 

0 

4 

1000 

1 

2 

Comer 

+ 

- 

- 

35 

4 

1000 

1 

3 

Comer 

- 

+ 

- 

0 

25 

1000 

3 

4 

Comer 

+ 

+ 

- 

35 

25 

1000 

3 

5 

Comer 

- 

+ 

0 

4 

17000 

1 

6 

Comer 

+ 

- 

+ 

35 

4 

17000 

1 

7 

Comer 

- 

+ 

+ 

0 

25 

17000 

3 

8 

Comer 

+ 

+ 

+ 

35 

25 

17000 

3 

9 

Face  Center 

- 

0 

0 

0 

14.5 

9000 

2 

10 

Face  Center 

+ 

0 

0 

35 

14.5 

9000 

2 

11 

Face  Center 

0 

- 

0 

17.5 

4 

9000 

1 

12 

Face  Center 

0 

+ 

0 

17.5 

25 

9000 

3 

13 

Face  Center 

0 

0 

- 

17.5 

14.5 

1000 

2 

14 

Face  Center 

0 

0 

+ 

17.5 

14.5 

17000 

2 

15 

Center 

0 

0 

0 

17.5 

14.5 

9000 

2 

16 

Center 

0 

0 

0 

17.5 

14.5 

9000 

2 

17 

0 

0 

0 

17.5 

14.5 

9000 

2 

18 

Center 

0 

0 

0 

17.5 

14.5 

9000 

2 

As  was  the  case  previously,  the  maximum  (+)  and  minimum  (-)  of  each  variable  are  used  as  well  as  the 
arithmetic  mean  of  these  values  (0)  to  design  the  experiments.  If  temperature  turns  out  to  be  important 
then  the  experiments  will  be  conducted  in  3  batches  that  correspond  to  &e  temperature  of  the  experiments 
in  order  to  save  time  and  reduce  costs.  This  will  eliminate  the  ability  to  block  or  randomize  Ae  results 
but  is  considered  a  necessary  step  given  the  large  impact  on  costs.  In  the  hypothetical  case  described 
here,  that  would  require  3  batches  of  5,  1 1,  and  5  experiments.  Analysis  of  four  center  point  replicates 
(runs  15-18)  will  provide  a  means  to  determine  the  pure  error  within  the  system. 

As  in  the  previous  experimental  designs  discussed,  the  dependent  variables  to  be  measured  in  the  course 
of  these  experiments  are  listed  in  Table  3-2.  Once  again,  only  those  dependent  variables  that  can  be 
quantified  can  be  analyzed  statistically.  That  is  to  say  the  particle  growth  r^  carmot  be  quantified. 

While  simple  ‘cube’  plots  or  X-Y  scatter  plots  can  provide  some  means  of  interpreting  the  results  of  a 
face-centered  cube  design,  generally  speaking  interpretation  will  require  computer  assistance.  We 
propose  to  use  IMP  Software  Version  3.1  fi'om  SAS  Institute,  Inc.  to  perform  our  analysis.  This  package 
allows  for  the  determination  of  the  b  coefficients,  standard  error  in  b,  and  t-values  for  the  quadratic 
solution  for  each  dependent  variable  being  considered.  The  results  can  be  graphically  depicted  using  2- 
dimensional  contour  plots. 
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In  the  case  of  our  hypothetical  3-&ctor  &ce-centered  cube  design  depicted  in  Table  3-7,  die  output  would 
include  calculation  of  each  of  the  b-coefBcients  for  each  factor  in  the  quadratic  formula; 


Y  =  bo+b,X14b2X2-rt)3X3+bi2XlX2+b,3XlX3+b23X2X3+buXl^-Hb22X2^-H)33X3^ 


Where; 


Y 

dependent  variable  of  interest  (Table  3-2) 

XI 

= 

temperature  (linear  effect) 

X2 

= 

salinity  (linear  effect) 

X3 

= 

bitumen  concentration  (linear  effect) 

X1X2 

=r 

temperature-salinity  (two-way  interaction  effect) 

X1X3 

= 

temperature-bitumen  concentration  (two-way  interaction  effect) 

X2X3 

= 

salinity-bitumen  concentration  (two-way  interaction  effect) 

Xl^ 

= 

temperature  (curvature  effect) 

X2^ 

= 

salinity  (curvature  effect) 

X3^ 

= 

bitumen  concentration  (curvature  effect) 

The  solution  to  this  equation  provides  a  means  of  predicting  the  outcome  for  any  “Y”  for  any  setting  of 
“X”.  The  fit  of  the  equations  for  each  “Y”  is  determined  by  standard  least  squares  regression  analysis. 

3.4  Utility  of  the  DOE  Phase  of  Experimentation 

The  outcome  of  the  benchtop  ejqjeriments  will  result  in 

(1)  identification  and  ranking  of  the  independent  variables  (Table  3-1)  which  have  the 
greatest  impact  on  each  of  the  dependent  variables  of  interest  (Table  3-2), 

(2)  a  narrowing/focussing  of  the  ranges  over  which  changes  in  the  independent  variables 
have  the  greatest  impact,  and 

(3)  qualitative  understanding  of  the  fete  of  Orimulsion  when  present  in  extremely  high 
concentrations  (50,000  mg/L). 

For  example,  if  salinity  is  determined  to  be  the  most  important  variable  in  determining  bitumen  particle 
size,  TPH  in  sediment,  or  dissolved  PAH,  our  benchtop  results  should  help  to  narrow  the  range  over 
which  salinity  is  most  critical,  e.g.,  15  to  20  7oo.  This  result  would  allow  us  to  design  a  focussed  set  of 
meso-scale  (flume)  experiments  to  address  particle  growth  over  this  narrow  range  of  salinity. 


4.  DESIGN  OF  BENCHTOP  EXPERIMENTS 

4.1  Overview  and  Description  of  Prototype  Benchtop  Apparatus 

As  described  in  Section  3.2,  the  fete  of  Orimulsion  spilled  onto  natural  waters  will  depend  on  the  complex 
interaction  of  physical,  chemical,  and  biological  processes.  Reproducing  the  effects  of  these  processes  in 
the  laboratory  is  extraordinarily  diflBcult,  essentially  impossible.  Keeping  this  in  mind,  the  benchtop 
experiments  to  be  performed  in  the  course  of  this  study  must  be  kept  relatively  simple  and  thereby 
allowing  for  some  reasonable  level  of  imderstanding  of  the  effects  of  fee  five  independent  variables 
(Table  3-1)  to  result.  The  design  of  fee  experiments  is  intended  to  provide  insight  into  the  expected 
behavior  of  Orimulsion  when  spilt  imder  a  wide  range  of  conditions.  The  objectives  of  fee  benchtop 
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Study  is  to  help  focus  the  meso-scale  experiments  that  will  follow  and  also  to  aid  in  the  design  of 
contingency  actions  under  a  wide  variety  of  conditions. 

In  this  section  we  describe  our  proposed  benchtop  design.  We  have  drawn  upon  our  first-hand  experience 
in  designing  and  testing  of  a  prototype  apparatus  (Fig.  4-1).  Our  design  has  attempted  to  incorporate 
practical  knowledge  obtained  in  previous  studies  on  Orimulsion  by  other  labs  (e.g.,  Intevep,  S.A.,  Univ. 
of  Mass.,  and  Univ.  of  Miami).  Some  key  elements  related  to  each  step  in  the  benchscale  experimental 
design  are  discussed  in  the  following  sections. 


Figure  4-1:  Photograph  showing  prototype  vessels  resting  on  and  seated  in  a  custom  made  orbital 
shaker  tabletop  capable  of  holding  four  vessels. 

4.2  Vessel  Construction 

Four-liter,  open  top,  clear  glass  vessels  will  be  fabricated  from  5  L  bottles  by  removing  their  tops  and 
polishing  the  cut  edges.  The  final  vessel  dimensions  will  be  24  cm  high  by  18  cm  diameter.  A  small 
sampling  port  consisting  of  a  glass  neck  containing  a  Teflon  stopcock  will  be  added  to  each  vessel.  This 
sampling  port  will  be  located  approximately  4  cm  from  the  bottom  of  the  vessels  (Fig.  4-1).  We  find  this 
location  for  the  sampling  port  provides  a  good  head  pressure  so  that  a  representative  water  sample  can  be 
drawn.  Samples  for  particle  size  analysis  (1,  3,  8,  and  24  hours),  TPH  in  water  (1,  3,  8  and  24  hours) 
will  be  collected  through  these  ports.  Attached  to  each  port  will  be  a  flexible  Tygon  tube  that  will  permit 
collection  of  these  samples  without  interruption  of  the  orbital  shaking  motion.  During  the  experiments 
the  vessels  will  be  loosely  covered  with  aluminum  foil  to  minimize  evaporation  and  contamination. 
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4.3  Salinity  Determination  and  Maintenance 

In  the  e?q)eiiments  to  be  conducted  the  salinity  of  die  water  will  be  determined  with  the  aid  of  a 
conductivity  meter.  Seawater  fixim  Duxbury  Bay  (~23  ppt)  will  be  adjusted  accordingly  using  Instant 
Ocean  sea  salts  or  freshwater  in  order  to  achieve  the  desired  salinity  for  each  experiment.  During  the  high 
temperature  experiments  (20-25°C)  there  may  be  some  change  in  the  salinity  of  the  water  due  to 
evaporative  losses.  However,  we  believe  that  this  change  will  be  minimized  by  die  use  of  aluminum  foU 
and  the  relatively  short  duration  of  these  experiments  (24  hours). 

4.4  Temperature 

The  benchtop  experiments  will  be  conducted  in  Battelle’s  temperature  controlled  rooms.  One  of  these 
rooms  contains  benches  for  the  benchtop  experiments  and  die  other  contains  the  4000  L  flume  tank.  Both 
rooms  must  be  operated  at  the  same  temperature.  The  range  of  temperatures  is  fiom  0°C  to  30°C. 

4.5  Orimulsion  Concentration 

A  new  shipment  of  Orimulsion  will  be  used  in  these  experiments  (see  below).  The  concentration  of 
Orimulsion  to  be  added  to  the  benchtop  experiments  will  be  determined  based  upon  the  measured  density 
of  this  material.  Widi  a  maximum  and  minimum  concentrations  of  Orimidsion  in  these  experiments 
being  17,000  mg/L  and  1000  mg/L,  respectively  (Table  3-1),  and  a  density  of  approximately  1.01  g/ml, 
we  can  expect  the  need  to  add  ^proximately  50.5  ml  and  2.97  ml,  respectively,  to  the  3  L  of  water  to 
obtain  the  desired  max.  and  min.  concentrations.'^  These  volumes  of  Orimulsion  will  be  precisely  added 
to  the  vessels  using  pre-weighed  50  ml  plastic  syringes.  Our  preliminary  experiments  have  shown  drat 
this  method  of  measuring  and  adding  Orimulsion  can  be  very  precise.  The  exact  volume  of  the 
Orimulsion  added  to  each  vessel  will  be  determined  by  weighing  the  syringes  after  squirting  out  as  much 
Orimulsion  as  possible  (as  some  will  stick  to  the  outside  of  the  syringe).  Orimrdsion  will  be  added  to  the 
vessels  after  the  shaker  tables  have  been  turned  on  in  order  to  assure  a  more  realistic  mixing  of  the 
product  in  the  water  column. 

4.6  Particulate  Load 

Piue  kaolinite  clay  has  been  obtained  from  the  Southeaster  Clay  Co.,  Aiken,  South  Carolina  (803-648- 
3246).  The  density  of  this  material  has  been  determined  to  be  2.6  g/ml.  The  particle  size  distribution  of  a 
kaolinite  slurry  has  been  measmed  using  the  Malvern  MasterSizer  X  (Section  5.2.2)  and  the  clay  was 
determined  to  have  a  bi-modal  particle  size  distribution.  The  first  population  was  centered  around  2 
microns  and  the  second  was  centered  around  25  microns.  It  is  our  impression  that  these  populations 
represent  the  thickness  (2  microns)  and  width  (25  microns)  of  these  platy  minerals.  The  total  particle  size 
ranged  ft^om  less  than  1.4  microns  to  100  microns.  The  source  and  characteristics  of  the  calcite  to  be  used 
has  not  yet  been  identified. 

The  mass  of  kaolinite  or  calcite  to  be  added  to  each  vessel  will  range  up  to  270  mg/L  which,  for  the  3  L 
benchtop  experiments,  indicates  810  mg  of  kaolinite  or  calcite  will  be  added.  We  propose  to  add  the 
kaolinite/calcite  to  the  vessels  and  vigorously  stir  the  water  to  achieve  a  complete  slurry.  Following  this 
the  orbital  shaker  table  will  be  turned  on  (see  Section  4.7)  and  imdoubtedly  some  of  the  clay  will  settle 
out,  certainly  imder  the  ‘mild’  agitation  conditions.  However,  some  will  remain  suspended  and  be  in 
suspension  at  the  time  when  the  Orimulsion  is  introduced.  We  believe  this  is  the  most  appropriate 
manner  by  which  to  treat  all  tire  experiments  equally.  The  initial  particle  size  distribution  of  the  clay 
slurry  will  be  determined  in  order  to  help  interpret  the  particle  size  results  once  the  Orimulsion  is  added  to 
the  system. 

The  four  extremely  hi^  concentration  experiments  (50,000  mg/L)  will  require  rpproximately  150  ml  Orimulsion 
be  added  tot  he  vessels  (Table  3-4’). 
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4.7  Energetics 

The  energetics  of  the  benchtop  ^paiatus  will  be  introduced  with  the  aid  of  a  modified  orbital  shaker  table 
(Figure  4-1).  Two  of  these  will  be  constructed  which  will  allow  for  8  experiments  to  be  performed  at  any 
given  time.  The  proposed  ‘mild’  and  ‘harsh’  agitation  conditions  are  10  and  100  rpm,  respectively.  Our 
preliminary  experiments  conducted  under  these  conditions  indicate  that  these  rpm  levels  should  provide 
an  adequate  range  of  turbulent  energy  in  the  3L  vessels.  Higher  rpm’s  tended  to  create  a  considerable 
vortex  within  the  vessel  which  acted  to  draw  Orimulsion  down  fiom  the  surfece  near  the  center  of  the 
vessel.  A  discussion  of  the  turbulent  energy  introduced  by  orbital  shaking  was  presented  in  Section 
3.2.I.4. 


5.  ANALYTICAL  METHODS  DURING  BENCHTOP  EXPERIMENTS 

5.1  Mixing  and  Characterization  of  Fresh  Orimulsion 

Prior  to  the  initiation  of  benchtop  testing,  the  new  shipment  of  Orimulsion  currently  en  route  to  Battelle 
will  be  thoroughly  mixed  using  an  air  powered  stainless  steel  drum  mixing  apparatus.  Aliquots  of 
approximately  4-L  will  be  drawn  off  into  one  gallon  paint  cans,  where  they  will  be  stored  for  fiiture  use 
throughout  these  benchtop  and  subsequent  meso-scale  flume  experiments.  The  new  “1998”  batch  of 
Orimulsion  will  be  put  through  an  abbreviated  suite  of  physical  and  chemical  measurements  to  ensure  the 
exact  makeup  of  the  material  relative  to  the  “1997”  Orimulsion  previously  analyzed  at  Battelle  (Ostazeski 
et  al.,  1998),  as  well  as,  any  specification  or  data  provided  by  Bitor.  These  measurements  will  include 
density,  viscosity,  latroscan  chemical  analysis  of  saturate,  aromatic,  resin,  and  asphaltene  (SARA) 
fractions,  and  bitumen  particle  size. 


5.1.1  Density 

The  density  of  the  flesh  Orimulsion®,  and  water-in-oil  emulsions  was  determined  using  an  Anton  Paar 
DMA  48  Density  Meter  oscillating  cell  digital  density  meter.  This  method  is  based  on  ASTM  D4052-91, 
Standard  Test  Method  for  Density  and  Relative  Density  of  Liquids  by  Digital  Density  Meter  (ASTM, 
1993).  Where  the  density  or  viscosity  of  the  bitumen  or  emulsions  formed  exceeded  the  practical 
operational  limits  of  the  density  meter,  the  density  will  be  obtained  gravimetrically  using  a  3  or  5  mL 
plastic  syringe.  The  mass  of  the  syringe  and  the  emulsion  is  then  made  and  the  mass  of  the  emulsion 
determined  by  difference  (equation  2.1). 

Mass  emulsion  (g)  =  [Mass  emulsion  (g)  +  Mass  syringe  (g)]  -Mass  syringe  (g) 

5.1.2  Viscosity 

The  viscosity  of  the  fresh  Orimulsion,  resulting  bitumen  and  any  water-in-oil  emulsions  will  be 
determined  using  a  Haake  RV-20  Rotational  Viscometer.  The  viscosities  will  be  determined  using  the 
viscometer,  operated  at  the  lowest  shear  rate  necessary  for  accurate  measurement  of  viscosity  (Daling  and 
Almas  1988).  Typically,  viscosities  for  flesh  and  weathered  oils  are  measured  at  shear  rates  of  1  s’’,  10  s’ 
',  100  s  ',  and  1000  s  '.  Based  on  our  past  experiences,  we  expect  the  viscosity  of  the  flesh  Orimulsion 
will  be  measured  at  100  s  '. 

5.1.3  SARA  Analysis 

The  relative  concentrations  of  saturated  hydrocarbons,  aromatic  hydrocarbons,  polars/resins  and 
asphaltenes  (SARA)  will  be  determined  by  latroscan  TLC/FID.  Approximately  10  pg  of  Orimulsion  or 
concentrated  water  extract  is  spotted  onto  a  silica  rod,  and  the  saturated,  aromatic,  and  resin  fractions  of 
the  oil  separated  by  development  of  the  rods  in  a  series  of  increasingly  polar  solvents  (w-hexane,  toluene. 
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dichloromethanermethanol).  The  relative  distribution  of  the  saturated,  aromatic,  and  resin  fractions  on 
each  rod  is  then  determined. 

5.1.4  Particle  Size  Analysis 

Particle  size  analysis  of  the  fresh  Orimulsion  will  conducted  using  a  Malvern  Mastersizer®  X  laser 
diffraction  particle  size  analyzer  (Fig.  5-1).  The  Mastersizer®  will  be  calibrated  initially  using  Duke 
Scientific  Corporation  (Palo  Alto,  CA)  NIST  particle  size  standards.  It  is  well  known  that  the  diameters 
of  non-spherical  particles  in  a  population  of  particles  can  be  calculated  and  reported  by  numerous 
methods.  These  methods  are  simply  different  ways  of  describing  the  diameters  of  non-spherical  particles, 
such  as  the  bitumen  particles  in  Orimulsion.  We  propose  to  use  the  most  commonly  referenced  type  of 
diameter,  the  volume  mean  diameter  (D[4,3]),  throughout  this  study.  The  volume  mean  diameter  is 
defined  as  the  diameter  of  a  sphere  having  the  same  volume  as  the  particles  being  measured. 

Reporting  a  single  diameter  for  a  broad  range  of  particles  can  be  misleading.  For  example,  if  multiple 
populations  of  particles  are  observed  the  volume  mean  diameter  for  the  entire  range  of  particles  is 
meaningless.  If  appropriate  the  raw  data  of  these  analyses  will  be  provided  so  that  the  particle  size 
distribution  is  available. 

Bitumen  particle  diameters  for  the  fresh  Orimulsion®  will  determined  by  adding  small,  well-stirred 
quantity  of  Orimulsion  to  a  circulating  volume  of  milli-Q  ultrapure  water  following  a  background 
measurements  on  the  water  alone.  An  obscuration  of  10-30%  will  be  obtained.  A  series  of  measurements 
will  be  obtained  and  averaged.  The  Malvern  MasterSizer  will  be  equipped  with  the  internal  flow  through 
sample  cell  (Fig.  5-1  A)  and  the  300  mm  focusing  lens.  This  lens  allows  for  measurement  of  particles  in 
the  1 .5  to  600  micron  range.  Processing  of  the  data  will  be  performed  using  the  polydispersion  model  and 
the  2$$D  presentation  code.  (Some  results  may  be  re-processed  and  reported  using  the  40HD 
presentation  code,  as  per  Intevep’s  method.  We  are  unable  to  use  a  1 18  mm  focal  lens  as  per  Intevep’s 
method). 


Figure  5-1:  Close-up  of  Mastersizer  X  internal  flow  through  cell  (A)  and  circulating  unit  (B). 

Orimulsion  will  be  added  to  circulating  (background)  Milli-Q  water  sample  through  the 
opening  in  the  circulating  unit  until  an  appropriate  obscuration  is  achieved.  Multiple 
measurements  will  immediately  follow  as  bitumen  particles  are  circulated  through  the  internal 
flow  through  cell,  thus  passing  through  the  laser  beam. 
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5.1.5  Gas  Chromatographic  Analysis  of  Fresh  Orimulsion 

Analysis  of  the  fiesh  Orimulsion*  for  Cg  to  C36  «-alkanes,  isoprenoid  hydrocarbons,  and  overall  boiling 
range  of  selected  samples  will  be  performed  according  to  modifications  of  EPA  Methods  8015.  The  total 
gas  chromatogr^hic  hydrocarbon  signature  of  fiesh  Orimulsion®  will  be  determined  using  a  Hewlett- 
Packard  5890  Series  n  c^illary  gas  chromatogrc^h  with  flame  ionization  detection  (GC/FID).  A  3  pL 
aliquot  of  a  5  mg/mL  dichloromethane  dilution  of  the  Orimulsion  will  be  injected  using  splitless 
techniques  onto  a  30-m  DB-5  cqiillary  column  (J&W  fused  silica  DB-5,  30-m,  0.32-mm  internal 
diameter,  and  0.25-|jm  film  thickness),  and  the  oven  will  be  prograrmned  fiom  35EC  to  320°C  at 
6°C/min.  Prior  to  sample  analysis,  a  multi-point  calibration  curve  containing  n-Cg  to  n-Cge  n-alkanes, 
pristane,  and  phjdane  will  be  generated  to  demonstrate  instrument  calibration  and  performance.  Based  on 
previous  experience  we  do  not  expect  any  resolved  alkanes  or  isoprenoids  to  be  present.  If  they  are, 
however,  they  will  be  quantified  by  the  method  of  internal  standards  using  the  internal  standard  5a- 
androstane. 

5.1.6  Polycyclic  Aromatic  Hydrocarbons  in  Fresh  Orimulsion 

The  analysis  of  target  parent  polycyclic  aromatic  hydrocarbons  (PAH’S),  alkylated  PAHs,  and 
dibenzothiophenes  (Table  5-1)  will  performed  on  the  fiesh  Orimulsion®  using  HRGC/MS  methods 
modified  fiom  EPA  Method  8270.  Basically,  prior  to  sample  analysis,  the  GC/MS  will  be  tuned  with 
perfluorotributylamine  (PFTBA),  and  a  multi-point  initial  calibration  consisting  of  parent  and  alkylated 
PAH  compounds  will  be  established  to  demonstrate  the  linear  range  of  the  analysis.  A  calibration  check 
standard  will  be  analyzed  in  order  to  monitor  instrument  response  relative  to  the  initial  calibration. 
Quantification  of  individual  components  will  be  performed  by  the  method  of  internal  standards  using 
relative  response  fiictors  (RRF)  for  individual  components  relative  to  the  internal  standards  acen^hthene- 
dio,  phenanthrene-dio,  and  benzo[a]pyrene-di2.  PAH  alkyl  homologues  will  be  quantified  using  the 
straight  baseline  integration  of  each  level  of  alkylation  and  the  RRF  of  the  respective  parent  PAH 
compound  (Table  5-1). 

The  fiesh  Orimulsion  sample  will  be  analyzed  by  GC/MS  using  a  Hewlett-Packard  5890  gas 
chromatograph  equipped  with  a  Hewlett-Packard  5970  mass  selective  detector  operating  in  the  selected 
ion  monitoring  (SIM)  mode.  A  2  pL  aliquot  of  the  sample  will  be  injected  into  a  gas  chromatogr^h 
equipped  with  a  high-resolution  capillary  column  (J&W  firsed  silica  DB-5-MS,  30-m,  0.25-mm  i.d.,  and 
0.25-pm  film  thickness)  operated  in  the  splitless  mode. 

5.1.7  Volatile  Organic  Compounds  (VOC) 

The  fresh  Orimulsion  will  also  be  analyzed  for  target  volatile  organic  compounds  (VOC;  Table  5-1)  using 
modified  EPA  8260  purge  and  trap  techniques,  utilizing  a  GC/MS  system  equipped  with  a  50-m 
Petrocol®  capillary  column.  Prior  to  sample  analysis,  the  instrument  will  be  tuned  with  PFTBA  and  a 
minimum  three-point  calibration  will  be  performed  to  demonstrate  the  linear  range  of  the  analysis.  The 
calibration  solution  is  composed  of  selected  VOCs  identified  in  Table  5-1. 
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Table  5-1.  Target  analyte  list  for  characterization  of  fresh  Orimulsion. 


GC/FID  Tareet  Analrtes 

SnUdns  Comnoimds  and  Renortins  Limits 

N^hthalene 

n-Cs  through  n-Cu 

GC/MS  SIS  Cottwounds 

Ci-naplidialenes 

Pristane 

Naphthalene-ds 

C2-napli1faaIettes 

Phytane 

Fluorene-dio 

Ca-napbdialenes 

Isoprenoid  1380 

Chiysene-di2 

C^-naphliialenes 

Isoptenoid  1470 

Phenol  -d5 

Biphenyl 

Isoprenoid  1650 

Acenaphthylene 

TPH 

GC/FID  SIS  Compound 

DibenzofiiTan 

o-teiphenyl 

Acenaphthene 

Fluorene 

Monoaromatic  Hvdrocarbons  (VOCsl 

GC/MS  SIS  Compounds 

Ci-fluotenes 

Benzene 

Acenaphthene-dio  ^ 

Qrfluorenes 

Toluene 

Phenanthrene-dio 

C3-fhiorenes 

Ethylbenzene 

Benzo[fl]pyrene-d)2 

Anthracene 

ni-Xylene 

Pfaenanthrene 

p-Xylene 

Ci'phenantfarenes/anthracenes 

o-Xylene 

Reportine  Luniis 

Cr-phenanthroies/anthracenes 

Isopropylbenzene 

PAH:  5  mg/kg 

C3-phenanthiienes/anthracenes 

»-Propylbenzene 

VOC:  10  mg/kg 

C4-phenanthienes/anthracenes 

l-Methyl-3-ethylbeozette 

alkanes;  100  mg/kg 

Dibenzothiophene 

l-Methyl-4-ethylbenzene 

TPH:  100,000  mg/kg 

Ci-dibenzothiophenes 

l-Methyl-2-ethylbenzene 

Cj-dibenzothiophenes 

1,3,5-Trimethylbenzene 

Cs-dibenzothiophenes 

1, 2,4/1, 2,3-Trimethylbenzene 

C4-dibenzothiophenes 

seC'Butylbenzene 

Fluoranthene 

l-Methyl-3-isopropylbenzene 

Pyrene 

1  -Methyl-4-isopropyIbenzene 

Ci-fluoranthenes/pyrenes 

l-Methyl-2-isopropyIbenzene 

C2-fluoranthenes/pyrenes 

l-Methyl-3-n-propylbenzene 

Cs-fluoranthenes/pyrenes 

1  -Methyl-4-n-propy  Ibenzene 

Benz[a]anthracene 

1 ,3,dimethyl-5-elhylbenzene 

Chrysene 

1 ,2-Diethylbenzene 

Cj-chrysenes 

1 ,3-Diethylbenzene 

Ci-chtysenes 

1,4-Diethylbenzene 

Cs-chiysenes 

l-Methyl-2-n-propylbenzene 

C4-chrysenes 

1 ,4-Dimethyl-2-ethylbenzene 

Benzo[b]fIuoranthene 

l,2-Diniethyl-4-ethylbenzene 

Benzo[A;]fIuoranthene 

1 ,3-Dimethy  1-2-ethylbenzene 

Benzo[e]pyrene 

l,3-Dimethyl-4-ethylbenzene 

Benzo[a]pyrene 

l,2-Dimethyl-3-ethylbenzene 

Peiylene 

1,2,4,5-Tetramethylbenzene 

lndeno[l,2,3-c,<flpyrene 

1,2,3,5-Tetramethylbenzene 

Dibenz[a,h]anthracene 

1,2,3,4-Tetramethylbenzene 

Benzo|;g,Ai]perylene 
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5.2  Collection  and  Analysis  of  Benchtop  Samples 

5.2.1  Collection  of  Samples  during  Benchtop  Experiments 

In  the  course  of  the  benchtop  experiments,  water  samples  will  be  collected  for  TPH  analysis  fiom  the 
vessel  sampling  ports  (Fig.  4-1)  at  intervals  of  1  and  8  hours.  The  orbital  shaker  tables  will  not  be  turned 
off  during  sampling  in  order  to  maintain  a  constant  degree  of  turbulence  in  the  water.  Approximately  40 
ml  water  samples  will  be  drawn  at  each  interval  from  each  sample  vessel  for  TPH  analysis.  Given  our 
experience  with  some  preliminary  experiments  the  reduced  quantity  of  water  collected  should  be 
sufficient  for  TPH  analysis  by  gravimetric  of  GC-FID  methods.  Some  difficulties  may  occur  during  the 
8-hour  sampling  when  under  some  circumstances  too  little  bitumen  might  be  still  suspended  in  the  water. 
In  these  cases  the  TPH  may  be  below  the  detection  limit.  The  water  sample  drawn  at  the  end  of  the 
experiments  (24  hr)  will  be  larger,  consisting  of  between  1  and  2  L.  This  volume  of  sample  will  permit 
the  sample  to  be  split  and  still  provide  enough  sample  to  obtain  good  concentration  data.  One  half  of  the 
water  will  be  used  for  TPH  and  total  PAH  analysis  and  one  half  will  be  used  for  dissolve  PAH  analysis. 
The  latter  aliquot  will  be  filtered  through  a  0.45  micron  filter  prior  to  analysis  in  order  to  remove  all  but 
the  dissolved  materials. 

At  the  1,  3,  8  and  24-hour  intervals,  samples  will  also  drawn  for  bitumen  particle  size  analysis.  It  should 
be  noted  (as  is  described  in  greater  detail  in  Section  5.2.2)  that  for  ffiose  experiments  that  include 
kaolinite  particulates,  it  will  not  be  possible  to  distinguish  bitumen  particle  size  firom  particulate  particle 
size.  In  the  latter  case,  particulates  will  include  bitumen,  clay,  and  bitumen/clay  conglomerates.  The 
volume  of  water  drawn  during  each  of  these  sampling  events  will  vary  slightly  depending  upon  the 
concentration  of  suspended  particulates  in  the  water.  Our  preliminary  experiments  have  shown  that  early 
in  the  experiments  as  little  as  10  ml  is  sufficient  to  achieve  appropriate  obscuration  levels  during  particle 
size  measurement.  However,  later  in  the  experiment,  when  very  little  suspended  particulate  is  present,  a 
larger  volume  100  ml  may  be  necessary.  Under  some  circumstances  we  can  imagine  that  not  enough 
suspended  matter  will  remain  in  the  water  to  achieve  a  reliable  particle  size  measurement.  In  such 
instances,  no  values  will  be  reported. 

5.2.2  Measurement  of  Particle  Size 

Measuring  particle  size  in  a  nuxture  can  be  done  by  a  variety  of  methods.  These  include  (1)  microscopy, 
(2)  sieving  and  weighing,  (3)  settling  experiments,  and  (4)  laser  diffraction.  We  believe  that  because  of 
the  potential  to  change  bitumen  particle  size  upon  changing  the  environment  of  the  sample,  only  laser 
diffraction  is  appropriate  for  the  bitumen/clay  particles  that  will  be  formed  in  the  course  of  these 
experiments.  However,  because  our  preliminary  experiments  have  indicted  that  particles  beyond  the 
range  of  the  Mastersizer  X  wiU  be  formed  (>  600  microns),  we  propose  to  supplement  the  laser  particle 
size  analysis  by  qualitative  microscopy. 

The  parameters  for  the  measurement  of  particulates  in  the  water  from  the  benchtop  experiments  using  the 
MasterSizer  X  will  be  the  same  as  described  for  fresh  Orimulsion  (Section  5.1.4).  Diameters  will  be 
reported  as  the  volume  mean  diameter  (D[4,3]),  which  refers  to  the  diameter  of  a  sphere  having  the  same 
volume  as  the  particles  being  measured.  Our  preliminary  experiments  have  indicated  that  there  may  be 
multiple  populations  present  in  any  given  sample.  If  this  is  the  case,  then  the  diameter  of  each  population 
will  be  determined  separately.  Qualitative  light  microscopy  will  be  used  to  describe  the  populations  (e.g., 
clay/calcite-dominant,  bitumen-dominant,  mixed,  etc.)  and  this  may  aid  in  determining  which  population 
is  most  representative  of  the  bitumen  that  is  present. 

As  noted  previously,  a  complication  in  the  measurement  of  bitumen  particle  size  lies  in  the  inability  to 
differentiate  bitumen  particles  from  clay/calcite  particles  in  the  samples.  This  will  be  problematic  as 
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particulate  load  is  considered  to  be  an  important  independent  variable.  We  believe  this  problem  can  be 
managed  as  the  real  goal  of  the  experiments  to  be  performed  is  to  recognize  changes  in  the  size  of  the 
bitumen  particles.  Our  preliminary  experiments  have  shown  that  even  with  kaolinite  in  the  system  the 
change  in  the  overall  particle  diameter  can  be  recognized  by  the  change  in  the  distribution  of  all  of  the 
particles  in  a  sample.  With  microscopy  to  supplement  the  laser  diffiaction  analysis  we  believe  a 
reasonable  interpretation  of  the  total  distributions  (i.e.,  bitumen  and  inorganic  particles)  can  be  made. 
However,  we  believe  that  it  would  be  inappropriate  to  use  the  average  diameters  of  these  mixed 
populations  in  any  attempt  to  ‘model’  bitumen  particle  growth  rate.  They  could  only  be  used  to  model  a 
total  particle  growth  rate.  (Of  course,  the  data  fiom  tiie  experiments  in  which  no  kaolinite  or  calcite  is 
added  will  provide  ‘bitumen-only’  particle  diameters,  and  these  results  could  be  used  to  model  a  bitumen 
particle  growth  rate). 

S.2.3  Extraction  of  Water  Samples  for  TPH  and  PAH 

Water  samples  collected  for  TPH  and  PAH  analysis  will  be  extracted  according  to  BOS  SOP  5-200.  Each 
batch  of  samples  will  include  a  procedural  blank  and  a  blank  spike.  Samples  will  be  spiked  with 
surrogate  internal  standards  (SIS)  prior  to  extraction  according  to  Table  5-2.  Final  extracts  will  be 
concentrated  to  a  volmne  of  0.5  mL.  All  extracts  will  then  be  spiked  with  recovery  internal  standard 
(RIS)  according  to  Table  5-3,  split,  and  transferred  to  tire  (jC/FID  and  GC/MS  task  leaders  for  TPH  and 
PAH  analysis,  respectively. 

Table  5-2.  Surrogate  Internal  Standard  (SIS)  spiking  information. 


Staridard  ID  ; 

Cbritents' 

Volume  to  Spike  (pL)  t 

j  EM82 

PAH/SHC  SIS 

50  i 

1  EM85  (BS  only) 

Alkanes 

50  I 

i  EM84  (BS  only) 

PAHs 

50  -  ! 

Table  5-3.  Recovery  Internal  Standard  (RIS)  spiking  information. 

Standard  CD " 

Contents 

Volume  to  Spike  (pL) 

1  EM83 

PAH/SHC  RIS 

!  50  1 

...  -  -  i 

5.2.4  Measurement  of  total  and  dissolved  TPH 

Analysis  of  water  sample  extracts  for  TPH  will  be  determined  by  HRGC  methods  modified  from  EPA 
Method  8015  and  defined  by  Battelle  SOP  5-202,  Determination  of  Low  Level  Total  Petroleum 
Hydrocarbon  and  Individual  Hydrocarbon  Concentrations  in  Environmental  Samples.  Prior  to  sample 
analysis,  a  minimum  of  a  five-point  calibration  will  be  performed  to  demonstrate  the  linear  range  of  the 
analysis.  The  calibration  solution  will  be  composed  of  selected  normal  and  iso-alkanes  ranging  in 
molecular  weight  fiom  Cg  through  C36.  Calibration  check  standards  will  be  analyzed  at  least  every  10 
samples  to  monitor  instnunent  response  relative  to  the  initial  calibration.  Instrument  calibration  and  other 
data  quality  objectives  are  detailed  in  SOP  202. 
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The  GC/FID  condition  for  the  analysis  will  be: 


Initial  Column  temperature: 

35°C 

Initial  hold  time; 

5  minutes 

Program  rate: 

6°C/mimrte 

Final  column  temperature: 

320°C 

Final  hold  time; 

10  minutes 

Injector  temperature: 

275°C 

Detector  temperature; 

325°C 

Column  flow  rate; 

1  mL/min 

The  reporting  limit  for  TPH  in  water  will  be  20.0  pg/L  (ppb). 

5.2.5  Measurement  of  total  and  dissolved  PAH 

The  analysis  of  the  target  PAH  compounds  will  be  performed  using  HRGC/MS  methods  modified  from 
EPA  Method  8270  and  defined  by  Battelle  SOP  5-157,  Identification  and  Quantitation  of  Polynuclear 
Aromatic  Hydrocarbons  (PAH)  by  Gas  Chromatography/Mass  Spectrometry.  Prior  to  sample  analysis, 
the  GC/MS  will  be  tuned  with  perfluorotributylamine  (PFTBA)  and  a  minimum  of  a  5-point  initial 
calibration  consisting  of  the  target  compounds  will  be  estabUshed  to  demonstrate  the  linear  range  of  the 
analysis. 

Analyte  concentrations  in  the  standard  solutions  will  range  fix>m  ^proximately  0.02  to  10  ng/pL. 
Quantification  of  individual  components  will  be  performed  by  the  method  of  internal  standards  using  the 
relative  response  factors  (RRFs)  of  the  parent  PAH  compounds  in  the  calibration  solutions.  Final  data 
will  be  reported  vs.  the  appropriate  surrogate  compound,  following  our  standard  IS  designation 
procedures  (only  the  PAH  data  will  be  surrogate  corrected).  Quantification  at  the  instrument  may  be  vs. 
the  RIS  acenaphthene-dio  and/or  fluorene-dio,  followed  by  surrogate  correction  of  the  data  in  a 
spreadsheet,  or  by  direct  quantification  vs.  the  surrogate  compound  in  the  instrument  method  whichever 
the  analyst  deems  most  efficient.  Instrument  calibration  and  other  data  quality  objectives  are  detailed  in 
SOP  5-157. 

The  GC/MS  condition  for  the  analysis  will  be: 


Initial  column  temperature:  40°C 

Irutial  hold  time  1  minute 

Program  rate:  6°/minutes 

Final  column  temperature:  290'’C 

Final  hold  time:  10  minutes 

Injector  temperature:  SOO^C 

Detector  temperature:  280°C 

Colunm  flow  rate:  1  mL/min 

Electronic  pressure  control  (EPC)  conditions: 

Vacuum  compensation:  On 

Pressure  at  injection:  15  psi 

Hold  time:  1 .50  min. 

Pressure  program  ramp:  99  psi/min. 

Final  pressure  7.7  psi 
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Reporting  limits  for  individual  PAHs  in  water  will  be  1.0  ng/L  (ppt). 

52.6  Measurements  for  Mass  Balance 

As  described  in  Section  3. 2.2.2,  considerable  knowledge  regarding  the  fate  of  Orimulsion  can  be  obtained 
by  determining  which  five  ‘compartments’  the  bitumen  ends  up  in  at  the  end  of  our  benchtop 
experiments: 

(1)  floating  bitumen/patty, 

(2)  sunk  bitumen/patty, 

(3)  dissolved  in  water, 

(4)  suspended  in  the  water,  or 

(5)  stuck  to  the  vessel  walls  (“bathtub  ring”). 

There  will  be  a  known  volume  of  Orimulsion  added  to  each  vessel  (determined  by  the  desired 
concentration).  Approximately  13%  of  this  volume  will  be  bitumen  and  the  density  of  bitumen  can  be 
used  to  determine  the  mass  of  bitumen  added  to  each  experiment.  Each  vessel  will  originally  contain  3  L 
of  water  and  immediately  following  the  addition  of  the  Orimulsion  all  of  the  bitumen  will  be  dispersed 
into  the  water.  After  24  hours  the  bitumen  will  occur  in  one  of  the  five  compartments. 

After  24  hours  any  floating  patty  (e.g..  Fig.  5-2)  will  be  recovered  using  pre-weighed  wooden  tongs  and 
placed  in  a  hood  to  dry  for  a  minimum  of  24  hours.  The  weight  of  the  floating  bitumen  will  then  be 
recorded. 


Figure  5-2:  Example  of  floating  surface  patty  (arrow)  formed  after  24  hours  in  a  preliminary 
experiment  involving  17,000  mg/L  Orimulsion  in  seawater.  Note  formation  of  the 
‘bathtub  ring’  of  bitumen  on  the  glass  walls  of  the  vessel. 
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After  24  hours  the  water  in  the  vessel,  a  slightly  reduced  volume  from  the  original  3L  due  to  the  1,  3  and 
8  hours  samplings,  will  be  drained  through  the  sampling  port  into  an  qjpropriate  container  in  preparation 
for  TPH  and  PAH  analysis.  The  volume  of  water  drained  will  be  recoMed.  The  dioroughly  mixed  water 
sample  will  be  split  jfor  dissolved  and  total  TPH  and  PAH  analysis  as  described  in  Sections  5.2.3  and 
5.2.4.  The  “dissolved”  split  will  be  filtered  through  pre-weighed  0.45  micron  filter  and  analyzed  to 
determine  the  concentration  of  dissolved  TPH  in  the  water  (Compartment  3).  The  mass  of  dissolved 
material  (Compartment  4)  will  be  calculated  given  the  total  volume  of  water  drained.  The  filter  will  be 
dried  and  weighed  in  order  to  determine  the  weight  of  total  (bitumen  and  perhrqps  clay)  suspended 
material  in  the  water.  The  “total”  water  split  will  be  extracted  and  analyzed  as  described  above.  The 
concentration  of  suspended  and  dissolved  TPH  will  be  determined  and  the  total  mass  will  be  calculated 
given  the  total  volume  of  water  drained.  The  mass  of  suspended  TPH  (Compartment  4)  will  be 
determined  by  the  difference  between  the  total  and  dissolved  masses. 

After  draining  the  water,  if  a  sinking  patty  is  present  it  will  be  collected,  dried,  and  weighed  as  described 
above.  If  no  patty  is  present,  only  bitumen  particles  these  will  be  recovered.  The  method  of  recovering 
this  may  vary  but  our  preliminary  e^qjeriments  have  shown  that  this  material,  mostly  a  mixture  of  bitumen 
and  clay,  can  be  rinsed  fiom  the  bottom  of  the  vessel  with  a  squirt  bottle  into  a  pre-weighed  beaker.  This 
bitumen/sediment  mixture  will  be  dried  and  weighed.  The  calculated  weight  of  the  sediment  (mineral 
matter)  will  be  subtracted  from  the  total  “sunk”  weight  to  yield  the  mass  of  the  sunk  bitumen 
(Compartment  2). 

Finally  the  now  “empty”  vessel,  having  been  pre-weighed  prior  to  the  experiment,  can  be  re-weighed  in 
order  to  determine  the  mass  of  bitumen  that  occurs  stuck  to  the  vessel  (Compartment  5).  Our  preliminary 
experiments  have  indicated  that  this  material  will  only  occur  at  the  water-air  internee  in  the  form  of  a 
‘bathtub  ring’. 

The  collective  total  of  the  masses  of  the  bitumen  in  each  compartment  should  ^proach  the  total  mass  of 
bitumen  spilled  in  each  experiment.  These  results  will  be  compared  for  the  given  experunental  conditions 
in  order  to  determine  the  influence  of  the  independent  variables  on  the  fate  of  Orimulsion. 

5.2.7  Reporting  of  Results 

All  analytical  data  generated  as  part  of  this  study  will  be  reviewed  by  Battelle’s  Quality  Assurance  Unit 
(QAU).  Data  will  be  subjected  to  a  targeted  audit  prior  to  delivery  to  the  client.  Data  will  be  reviewed  by 
the  Project  Managers  for  chemical  reasonableness.  All  project  deliverables  will  be  subjected  to  peer 
review  for  technical  and  editorial  content.  Prior  to  issuance  toe  final  Phase  11  and  Phase  III  reports  and 
supporting  documentation  (e.g.  relevant  figures,  tables  and  data)  will  undergo  a  QA  audit. 

A  report  siunmarizing  toe  results  of  toe  benchtop  experiments  will  be  issued  at  toe  completion  of  Phases 
II  and  III.  The  Phase  II  report  will  include  a  discussion  of  toe  effects  of  toe  five  independent  variables 
(Table  3-1)  on  toe  various  dependent  variables  (Table  3-2)  and  on  toe  mass  balance  results  as  determined 
from  toe  benchtop  studies.  This  report  will  also  include  toe  results  of  our  flume  energetic  study  and  our 
recommendations  for  Phase  Hi’s  meso-scale  flume  experiments  generally  described  below. 


6.  PHASE  III  -  MESO-SCALE  FLUME  TANK  EXPERIMENTS 

6.1  Introduction 

The  meso-scale  experiments  comprising  Phase  III  of  this  study  will  be  conducted  in  Battelle’s  4000  L 
flume  tank  under  appropriate  energy  conditions  (see  below).  The  specifications  for  toe  tank  have  been 
described  elsewhere  (Fredriksson  et  al.,  1996).  The  tank  will  be  prepared  for  toe  weathering  experiments 
by  filling  it  with  toe  desired  amount  and  type  of  water.  The  initial  water  to  be  used  will  be  (1)  seawater 
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from  Duxbury  Bay,  MA  that  will  be  acquired  with  Battelle’s  seawater  system  arid/or  (2)  fresh  water 
obtained  from  Battelle’s  fresh  water  well.  The  tar^t  salinity  will  be  obtained  by  mixing  these  two  water 
sources  and/or  by  the  addition  of  Instant  Ocean™  sea  saJts.  The  salinity  will  be  monitored  with  a 
conductivity  meter  calibrated  for  the  temperature  of  the  water.  Both  seawater  and  freshwater  will  be 
filtered  to  5  microns  prior  to  its  use  in  the  flume  experiments. 

The  temperature  conditions  that  are  ultimately  decided  for  the  flume  experiments  will  be  set  for  the  flume 
room.  The  temperature  of  the  water  added  to  the  tank  will  be  monitored  until  it  reaches  the  temperature 
set  for  the  room. 


Figure  6-1:  Battelle’s  4000  L  flume  tank  showing  conditions  producing  spilling  breaking  waves. 

The  flume  has  a  track  length  of  9.1  m  and  a  depth  of  0.8  m.  Fifteen  wind  vanes  can  generate 
up  to  25  m/sec  wind  speed  and  the  wedge-shaped  plunger  (foreground)  stroke  ranges  from  2- 
30  cm  at  rates  up  to  100  cpm. 

6.2  Sample  Collection  and  Preparation 

Once  the  flume  experiments  are  initiated,  samples  of  Orimulsion  within  the  water  column  will  be  taken  at 
pre-set  times  for  analysis  according  to  the  schedule  and  analyses  summarized  in  Table  6-1 .  Naturally,  the 
schedule  and  analysis  may  be  changed  after  considering  the  Phase  II  results.  The  parameters  to  be 
measured  in  the  course  of  these  experiments  will  likely  parallel  those  made  at  the  benchtop  scale  (Table 
3-2)  and  those  from  previous  flume  experiments  on  Orimulsion  (Ostazeski  et  al.,  1998). 

At  the  indicated  sampling  times  following  introduction  of  the  Orimulsion,  water  samples  will  be  collected 
using  a  stainless  steel  hand  pump.  The  samples  will  be  collected  from  a  depth  of  approximately  35  cm 
below  the  water  surface  and  mid  channel  in  the  flume  track.  The  water  column  samples  will  be  prepared 
for  the  appropriate  analyses  according  to  standard  operating  procedures  described  in  Section  5-2. 
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Table  6-1.  Possible  Sampling  and  Testing  Scheme  for  Orimulsion'^  during  Phase  III  Flume  Testing 


Fhtatmg^Sinldng  IKtumen  Unfiltered  Water  Filtered  Water 


Time 

(hours) 

Viscosity 

Density 

TPH 

Analysis 

(sediment) 

SARA 

Analysis 

(latroscan) 

TPH 

Analysis 

PAH 

Analysis 

Particle 

Size 

TPH 

Analysis 

PAH 

Analysis 

0 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

15  min 

— 

— 

_ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

30  min 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

1 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

2 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

4 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

8 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

12 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

24 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

48 

— 

_ 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

72 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

96 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

120 

— 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

144 

-  . 

— 

— 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

168 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

6.3  Phase  nC  -  Evaluation  and  Quantification  of  Energy  within  the  Flume  Tank 

The  quantification  of  the  energy  conditions  prior  to  the  spilling  of  any  Orimulsion  in  the  flume  tank  is  a 
separate  and  integral  part  of  this  study.  This  will  be  particularly  true  if  the  benchtop  studies  reveal  that 
energetics  plays  a  very  important  role  in  the  fete  of  Orimulsion  spilled  on  water.  The  problem  remains 
separate,  however,  because  translating  the  kinetic  energy  levels  fi-om  the  benchtop  to  the  meso-scale  will 
be  essentially  impossible. 

The  conditions  of  the  earlier  Orimulsion  studies  (Ostazeski  et  al.,  1997a;  1997b;  1998)  were  arbitrarily  set 
to  match  those  previously  used  by  IKU  (Singsaas  et  al.,  1995).  The  ‘energetics’  of  the  previous  studies 
conducted  in  Battelle’s  flume  taitk  were  all  performed  under  these  conditions  which  were  described  by 
Fredriksson  et  al.  (1996).  In  that  study  the  following  data  were  generated;  (1)  wave  height  measurements 
using  a  c^acitance  wave  probe  and  (2)  current  velocity  measurements  at  the  center  of  the  channel  at  six 
depths  using  acoustic  Doppler  velocimeters  (ADV)  and  neutrally-buoyant  styrene  beads.  These 
measurements  were  made  with  the  plunger  and  wind  vane  operating  at  the  same  levels  as  in  the  previous 
Orimulsion  spill  experiments,  i.e.,  wind  vanes  at  80  cpm  and  plunger  at  50  cpm  with  a  12  cm  stroke 
(Ostazeski  et  al.,  1997a;  1997b;  1998).  The  energy  density  under  these  conditions  was  determined  to  be 
22.39  Joules/m^,  which  is  well  within  the  energy  density  estimates  for  open  ocean  sea-state  2  conditions 
(10.48  to  298.61  Joules/m^).*  Stokes  second  order  wave  theories  were  used  to  determine  the  mean 
vertical,  radial  and  tangential  current,  wave  and  turbulent  velocity  components  imder  these  conditions. 

Concern  was  recently  expressed  that  the  earlier  work  performed  on  Orimulsion,  in  feet,  on  all  previously 
spilled  oils,  were  conducted  at  energies  that  were  too  high.  The  basis  for  this  concern  was  the  feet  that  in 
the  event  of  a  real  spill,  there  is  normally  a  decrease  in  the  turbulent  mixing  with  depth  (P.  Masciangioli, 
1997,  personal  communication).  However,  it  was  observed  that  in  the  flmne  comparable  turbulence 

'  Note  that  without  the  plunger  operating  the  energy  was  only  7.06  Joules/m^. 
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existed  throughout  the  entire  flume  water  depth  leading  to  the  complete  mixing  of  the  water  mass.  This 
was  suggested  to  be  the  reason  that,  during  the  fieshwater  experiments,  the  bitumen  was  not  observed  to 
settle  out  of  die  water.  Perh^s  the  simplest  solution  to  this  concern  would  be  to  eliminate  the  use  of  the 
plunger  ^paratus  since  it  is  presumably  responsible  for  introducing  most  of  the  energy  into  the  deeper 
parts  of  the  flume  tank. 

The  “flume  energetics”  aspect  of  the  proposed  study  undoubtedly  reqmres  additional  discussions  between 
Bitor,  NOAA,  USCG  and  Battelle  persormel.  We  anticipate  the  scope  of  this  phase  of  the  study  to  be 
determined  collectively  during  our  meeting  March  19,  1998  in  Duxbury.  Both  the  technical  and  financial 
limitations  of  the  study  must  be  kept  in  mind  throughout  these  discussions. 

It  is  our  original  proposal  to  perform  a  series  of  measurements  in  the  flume  tank  under  a  variety  of  wind 
vane  and  plunger  settings.  We  propose  to  measure; 

(1)  wave  height, 

(2)  wave  period,  and 

(3)  wind  speed. 

It  is  our  opinion  that  these  gross  measurements  may  allow  for  an  easier  ‘translation’  of  real  world 
conditions  and  if  particle  velocities  were  measured  and  modeled.  Therefore,  we  do  not  intend  to  measure 
current  velocity  nor  calculate  current,  wave  or  turbulence  component  velocities.  Instead,  we  would 
simply  use  the  matrix  of  wind  vane/plimger  conditions  versus  wave  height,  wave  period  and  wind  sped 
data  to  select  the  appropriate  wind  vane  and  plunger  settings  to  be  used  for  the  meso-scale  studies.  We 
would  propose  to  acquire  this  data  at  the  same  time  that  the  benchtop  studies  would  be  being  conducted. 
We  would  seek  agreement  from  Bitor  and  USGC  representative  on  the  conditions  imder  which  the  meso- 
scale  studies  should  be  conducted  in  order  to  more  closely  match  real  world  conditions. 

After  receiving  additional  comments  on  this  workplan,  we  have  agreed  to  consider  including  a  study  by 
which  the  turbulent  dififiisity  through  the  flume  can  be  investigated  for  a  small  munber  of  wind/wave 
conditions  or  settings.  The  data  from  these  investigations  could  be  used  by  NOAA  and  Intevep 
representatives  to  aid  in  determining  the  wind/wave  conditions  under  which  to  conduct  the  flume 
experiments  and  allow  for  modeling  the  results  in  open  water  conditions.  This  study  will  involve  the  use 
of  inert  particles  of  a  known  diameter  and  density.  (The  nature  of  these  materials  is  yet  to  be  determined 
and  we  expect  to  receive  some  guidance  from  NOAA  on  an  appropriate  product  to  use)  Let’s  assume 
these  are  slightly  buoyant  particles.  The  buoyant  particles  will  be  added  to  the  flume  tank  for  a  given 
wind/wave  condition.  The  turbulent  energy  within  the  tank  will  drive  them  down  into  the  tank  and 
distribute  the  particles  throughout  the  water  colutim.  Buoyant  forces  will  attempt  to  drive  them  back  up, 
but  turbulence  will  keep  some  down.  These  forces  should  result  in  a  fairly  “steady  state”  at  which  the 
concentration  of  particles  decreases  with  increasing  depth.  Using  video  or  photogr^hy,  the  number  of 
particles  at  a  given  depth  range  will  be  determined  as  accurately  as  possible.  The  number  density  of 
particles  at  discrete  depths  within  the  tank  can  be  used  to  calculate  the  turbulent  difiiisivity  according  to 
the  following  equation; 


dN/dt  =  D  (d^N/dz^)  -  U  (dN/dz) 

Where  D=turbulent  diffusivity,  N  =  number  density  of  particles,  n  =  number  of  particles,  U  =  Stokes 
velocity,  z  =  depth. 

If  still  considered  to  be  important  enough  to  justify  the  additional  costs,  and  some  technical  guidance 
becomes  available  regarding  the  selection  of  ‘inert  particles’,  we  will  attempt  to  collect  number  density 
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for  a  small  number  of  wind/wave  conditions  that  will  allow  some  cqiacity  to  translate  the  flume 
eneigetics  to  open  water  energetics.  Battelle’s  role  in  this  phase  of  the  study  is  limited  to  collection  of 
data  imder  guidance  from  NOAAAJSGC  and  Bitor/Intevep  representatives.  The  costs  for  this  additional 
component  of  the  study  are  not  included  in  Section  7. 


7.  TIME  CONSIDERATIONS 

This  workplan  constitutes  the  deliverable  for  Phase  I  of  this  project.  Table  7-1  summarizes  the  tasks  that 
will  be  performed  in  both  Phases  II  and  m.  Because  of  the  dynamic  character  of  Phases  n  and  III 
Battelle  wiU  provide  monthly  updates  regarding  the  progress  and  costs  accmed  throughout  the  project. 
This  will  allow  constant  re-assessment  of  the  direction  and  scope  of  subsequent  phases. 

Table  7-1:  Inventory  of  Tasks  to  be  Performed  in  each  Phase  of  the  Study 


PHASES  AND  TASKS 

Phase  I  -  Design  of  Behavior  Characterization  Experiments 

Planning  Meeting 

Literature  Review 

Technical  Meetings 

Evaluation  of  particle  size  analyzers 

Development  of  Design  of  Experiment 

Development  Benchtop  Prototype 

Conducting  Preliminary  Benchtop  Experiments 

Preparation  of  Phase  I  deliverable 

^  ■ 

Phase  II  -  Evaluate  Fate  and  Determine  Conditions  for  Meso-Scale  Testing 

Phase  nA  -  Conduct  Benchtop  Experiments 

analytical  costs  -  fresh  Orimulsion 

constract  benchtop  apparatus,  purchase  materials 

perform  12  screening  design  experiments 

Report  screening  results 

perform  up  to  47  three-level  response  surface  design  experiments 

analytical  costs  -  dependent  variables  and  mass  balance 

Phase  ns  -  Data  Interpretation 

DOE  evaluation  using  IMP  software 

Phase  UC  -  Energetic  Study 

Constmct  wind  speed/wave  height/wave  frequency  matrix 

Phase  11  Report  Preparation  with  Recommendations  for  Phase  III 

Review  Meeting  -  Duxbury 

Phase  III  -  Meso-Scale  Flume  Tank  Tests 

Conduct  4  or  5  meso-scale  flume  tests 

Analytical  Costs  -  dependent  variables  of  interest  in  Phase  11 

Phase  ni  Report  with  Final  Section  on  Synthesis  of  Project 
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APPENDIX  III 


Response  Surface  Design  Model 
Calctilations 


»creenlng  Fit 
(PS1  ) 


(Summary  of  Fit ) 

' 

RSquare 

0.759198 

RSquare  Adj 

0.548497 

Root  Mean  Square  Error 

42.63504 

Mean  of  Response 

73.49548 

Observations  (or  Sum  Wgts) 

31 

(Analysis  of  Variance"^ 


Source 

DF 

Sum  of  Squares 

Mean  Square 

Model 

14 

91695.78 

6549.70 

Error 

16 

29083.95 

1817.75 

C  Total 

30 

120779.73 

(Lack  of  Fit  } 
(Parameter  Estimated) 


Term 

Estimate 

Std  Error 

t  Ratio 

Intercept 

40.498531 

41.96349 

0.97 

Salinity 

6.2048185 

3.297419 

1.88 

OR400 

0.0023283 

0.007949 

0.29 

Energy 

-2.864445 

1.516547 

-1.89 

Loading 

1.297812 

0.427443 

3.04 

Salinity'Salinity 

-0.170765 

0.086419 

-1.98 

OR400*Salinity 

0.0000758 

0.000076 

1.00 

OR400*OR400 

-1e-7 

4.135e-7 

-0.24 

Energy*Salinity 

0.0073556 

0.013535 

0.54 

Energy*OR400 

0.0000081 

0.00003 

0.27 

Energy*Energy 

0.0238164 

0.01307 

1.82 

Loading*Salinity 

-0.001951 

0.004512 

-0.43 

LQading*OR400 

-0.000003 

0.00001 

-0.30 

Loading*Energy 

0.0003928 

0.001755 

0.22 

Loading*Loading 

-0.004774 

0.001452 

-3.29 

Prob>|tl 

0.3489 

0.0782 

0.7734 

0.0772 

0.0079 

0.0657 

0.3340 

0.8119 

0.5943 

0.7877 

0.0872 

0.6711 

0.7683 

0.8257 

0.0046 


Greening  Fit 


;PS3  1 


(Summary  of  Fit ) 

•  - 

RSquare 

0.743712 

RSquare  Adj 

0.519461 

Root  Mean  Square  Error 

46.06573 

Mean  of  Response 

65.05903 

Observations  (or  Sum  Wgts) 

31 

i(Analysis  of  VariancQ 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Model 

14 

98526.52 

7037.61 

3.3164 

Error 

16 

33952.83 

2122.05 

Prob>F 

C  Total 

30 

132479.34 

0.0121 

V. 


[(Lack  of 

(Parameter  Estimater^ 


Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

48.721588 

45.34014 

1.07 

0.2985 

Salinity 

12.502851 

3.562751 

3.51 

0.0029 

OR400 

0.0170416 

0.008589 

1.98 

0.0647 

Energy 

-4.833287 

1.638578 

-2.95 

0.0094 

Loading 

0.2018765 

0.461838 

0.44 

0.6679 

Saiinity*Salinity 

-0.353135 

0.093373 

-3.78 

0.0016 

OR400*Salinity 

0.0000479 

0.000082 

0.58 

0.5682 

OR400‘OR400 

-9.132e-7 

4.468e-7 

-2.04 

0.0578 

Energy*Salinity 

-0.000777 

0.014624 

-0.05 

0.9583 

Energy*OR400 

0.0000005 

0.000032 

0.02 

0.9869 

Energy*Energy 

0.041263 

0.014121 

2.92 

0.0100 

Loading*Salinity 

-0.001479 

0.004875 

-0.30 

0.7656 

Loading*OR400 

-0.000003 

0.00001 1 

-0.29 

0.7770 

Loading*Energy 

0.0007787 

0.001896 

0.41 

0.6867 

1  Loading*Loading 

-0.000603 

0.001569 

-0.38 

0.7059 

(Parameter  Estimates 

D 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

32.095051 

51.08847 

0.63 

0.5387 

Salinity 

12.773623 

4.014445 

3.18 

0.0058 

OR400 

0.0263109 

0.009677 

2.72 

0.0152 

Energy 

-3.596207 

1.846321 

-1.95 

0.0692 

Loading 

-0.58055 

0.520391 

-1.12 

0.281 1 

Salinity'Salinity 

-0.348884 

0.105211 

-3.32 

0.0044 

OR400*Salinity 

0.0000195 

0.000093 

0.21 

0.8364 

OR400*OR400 

-0.000001 

5.035e-7 

-2.75 

0.0142 

Energy*Salinity 

-0.002492 

0.016478 

-0.15 

0.8817 

Energy*OR400 

-0.00001 1 

0.000036 

-0.30 

0.7650 

Energy*Energy 

0.0312045 

0.015912 

1.96 

0.0675 

Loading*Salinity 

-0.001763 

0.005493 

-0.32 

0.7523 

Loading*OR400 

-0.000003 

0.000012 

-0.25 

0.8036 

Loading*Energy 

0.0006202 

0.002136 

0.29 

0.7753 

Loading*Loading 

L 

0.0019681 

0.001768 

1.11 

0.2821 

Effect  Test 


Prediction  Profile 


(Parameter  Estimates  ) 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

26.35642 

66.24229 

0.40 

0.6960 

Salinity 

10.659205 

5.205206 

2.05 

0.0574 

OR400 

0.0103664 

0.012548 

0.83 

0.4209 

Energy 

-4.831103 

2.393975 

-2.02 

0.0607 

Loading 

0.9624622 

0.674749 

1.43 

0.1730 

Salinity*Salinity 

-0.277094 

0.136418 

-2.03 

0.0592 

OR400*Salinity 

-0.00003 

0.00012 

-0.25 

0.8073 

OR400*OR400 

-4.154e-7 

6.528e-7 

-0.64 

0.5336 

Energy*SaIinity 

-0.000165 

0.021366 

o 

b 

0.9939 

Energy*OR400 

-0.00003 

0.000047 

-0.65 

0.5245 

Energy'Energy 

0.0462938 

0.020631 

2.24 

0.0393 

Loading*Salinity 

-0.001607 

0.007122 

-0.23 

0.8244 

Loading*OR400 

-0.000005 

0.000016 

-0.32 

0.7512 

Loading*Energy 

0.0014041 

0.00277 

0.51 

0.6191 

Loading*Loading 

-0.003545 

0.002292 

-1.55 

0.1416 

Screening  Fit 


Prediction  Profile 


Salinity 


OR400 


Energy 


Loading 


270 


icreening  Fit 


. . . 

(Summary  of  Fit ) 

RSquare 

0.793133 

RSquare  Adj 

0.612124 

Root  Mean  Square  Error 

1.805823 

Mean  of  Response 

2.874194 

Observations  (or  Sum  Wgts) 

31 

(Analysis  of  Variance^ 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Model 

14 

200.04338 

14.2888 

4.3817 

Error 

16 

52.17597 

3.2610 

Prob>F 

C  Total 

30 

252.21935 

0.0030 

V  _ / 

[(Lack  of  nf) 


(Parameter  Estimates^ 


Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

0.0374578 

1.777379 

0.02 

0.9834 

Salinity 

0.0861456 

0.139663 

0.62 

0.5460 

OR400 

0.0001886 

0.000337 

0.56 

0.5831 

Energy 

-0.024585 

0.064234 

-0.38 

0.7069 

Loading 

0.0119271 

0.018105 

0.66 

0.5194 

Salinity*Salinity 

-0.00361 

0.00366 

-0.99 

0.3387 

OR400*Salinity 

0.0000089 

0.000003 

2.75 

0.0141 

OR400*OR400 

-7.898e-9 

1.752e-8 

-0.45 

0.6581 

Energy*Salinity 

0.0011349 

0.000573 

1.98 

0.0652 

Energy*OR400 

0.0000026 

0.000001 

2.04 

0.0588. 

Energy*Energy 

0.0001701 

0.000554 

0.31 

0.7625 

Loading*Salinity 

-0.000294 

0.000191 

-1.54 

0.1439 

Loading*OR400 

-7.81 3e-7 

4.18e-7 

-1.87 

0.0800 

Loading*Energy 

-0.000135 

0.000074 

-1.81 

0.0885 

Loading*Loading 

0.0000079 

0.000062 

0.13 

0.8990 

V _ / 


creening  Fit 


(Summary  of  Fit ) 

RSquare 

0.913948 

RSquare  Adj 

0.838653 

Root  Mean  Square  Error 

0.127845 

Mean  of  Response 

0.384694 

Observations  (or  Sum  Wgts) 

31 

(Analysis  of  Variance  J 


Source 

DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

Model 

14 

2.7774886 

0.198392 

12.1382 

Error 

16 

0.2615111 

0.016344 

Prob>F 

iC  Total 

30 

3.0389997 

<.0001 

(Lack  of  Fit  ] 


(Parameter  Estimates  ) 

< 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

Intercept 

-0.155355 

0.125832 

-1.23 

0.2348 

Salinity 

0.033053 

0.009888 

3.34 

0.0041 

OR400 

0.0000322 

0.000024 

1.35 

0.1953 

Energy 

0.0131156 

0.004548 

2.88 

0.0108 

Loading 

-0.002399 

0.001282 

-1.87 

0.0796 

Salinity*Salinity 

-0.000335 

0.000259 

-1.29 

0.2147 

OR400*Salinity 

0.0000004 

2.283e-7 

1.70 

0.1082 

OR400‘OR400 

-2.038e-9 

1.24e-9 

-1.64 

0.1197 

Energy*Salinity 

-0.000173 

0.000041 

-4.27 

0.0006 

Energy*OR400 

0.0000002 

8.878e-8 

2.48 

0.0247 

Energy*Energy 

-0.000126 

0.000039 

-3.21 

0.0055 

Loading*Salinity 

-0.000007 

0.000014 

-0.49 

0.6341 

Loading*OR400 

-2.764e-8 

2.959e-8 

-0.93 

0.3641 

Loading^Energy 

-0.000005 

0.000005 

-0.93 

0.3653 

Loading*Loading 

0.00001 

0.000004 

2.29 

0.0359 

_ : _ ; 


Greening  Fit 


(Summary  of  Fit 


fAnalysis  of  Variance 
Source  DF  Su 

Model  14 

Error  16 

C  Total  30 


DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

14 

2.6477745 

0.189127 

8.4624 

16 

0.3575858 

0.022349 

Prob>F 

30 

3.0053603 

<.0001 

CParameter  Estimates  ) 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>ll| 

Intercept 

0.5019611 

0.147141 

3.41 

0.0036 

Salinity 

-0.03787 

0.011562 

-3.28 

0.0048 

OR400 

-0.000051 

0.000028 

-1.84 

0.0848 

Energy 

0.0064977 

0.005318 

1.22 

0.2394 

Loading 

0.0020299 

0.001499 

1.35 

0.1944 

Salinity*Salinity 

0.0006596 

0.000303 

2.18 

0.0448 

OR46o*Salinity 

-6.663e-7 

2.67e-7 

-2.50 

0.0239 

OR400*OR400 

4.241 5e-9 

1.45e-9 

2.93 

0.0099 

Energy*Salinity 

0.0001224 

0.000047 

2.58 

0.0202 

Energy*OR400 

-1.057e-7 

1.038e-7 

-1.02 

0.3237 

Energy*Energy 

-0.000097 

0.000046 

-2.11 

0.0505 

Loading*Salinity 

-0.000004 

0.000016 

-0.26 

0.7971 

Loading*OR400 

-6.907e-8 

3.461  e-8 

-2.00 

0.0633 

Loading*Energy 

0.000019 

0.000006 

3.09 

0.0070 

Loading*Loading 

-0.000005 

0.000005 

-1.00 

0.3335 

tilllllll 


icreening  Fit 


C%Ring  ] 


[Summary  of  Fit  ] 

RSquare 

0.855077 

RSquare  Adj 

0.728269 

Root  Mean  Square  Error 

0.107847 

Mean  of  Response 

0.268977 

Observations  (or  Sum  Wgts) 

31 

[Analysis  of  Variance 
Source  DF  Su 

Model  14 

Error  16 

C  Total  30 


DF 

Sum  of  Squares 

Mean  Square 

F  Ratio 

14 

1,0980061 

0.078429 

6.7431 

16 

0.1860966 

0.011631 

Prob>F 

30 

1.2841027 

0.0003 

[Parameter  Estimates  ) 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>|tl 

Intercept 

0.4306283 

0.106149 

4.06 

0.0009 

Salinity 

-0.002011 

0.008341 

-0.24 

0.8125 

OR400 

-0.000001 

0.00002 

-0.05 

0.9596 

Energy 

-0.004101 

0.003836 

-1.07 

0.3009 

Loading 

-0.000951 

0.001081 

-0.88 

0.3921 

Salinity*Salinity 

-0.000043 

0.000219 

-0.20 

0.8463 

OR400*Salinity 

-5.116e-8 

1.926e-7 

-0.27 

0.7939 

OR400*OR400 

-7.73e-10 

1.046e-9 

-0.74 

0.4709 

Energy*Salinity 

0.0000437 

0.000034 

1.28 

0.2196 

Energy*OR400 

-1.282e-7 

7.489e-8 

-1.71 

0.1062 

Energy*Energy 

0.0000853 

0.000033 

2.58 

0.0201 

Loading*Salinity 

0.0000062 

0.00001 1 

0.54 

0.5956 

Loading*OR400 

8.6285e-8 

2.496e-8 

3.46 

0.0033 

Loading'Energy 

-0.000018 

0.000004 

-3.95 

0.0011 

Loading*Loading 

0.0000011 

0.000004 

0.29 

0.7745 

Effect  Test 


Prediction  Profile 


^Parameter  Estimates 

Term 

Intercept 

Salinity 

OR400 

Energy 

Loading 

Salinity*Salinity 

OR400*Salinity 

OR400*OR400 

Energy'Salinity 

Energy*OR400 

Energy*Energy 

Loading*Salinity 

Loading*OR400 

Loading*Energy 

Loading*Loading 


Estimate 

Std  Error 

t  Ratio 

Prob>|t| 

0.2218206 

0.079116 

2.80 

0.0127 

0.0068229 

0.006217 

1.10 

0.2887 

0.0000202 

0.000015 

1.35 

0.1961 

-0.015513 

0.002859 

-5.43 

<.0001 

0.0013159 

0.000806 

1.63 

0.1220 

-0.000281 

0.000163 

-1.73 

0.1036 

0.0000003 

1.435e-7 

2.29 

0.0362 

-1.439e-9 

7.8e-10 

-1.85 

0.0836 

0.0000069 

0.000026 

0.27 

0.7906 

1.3576e-8 

5.582e-8 

0.24 

0.8109 

0.0001374 

0.000025 

5.58 

<.0001 

0.0000045 

0.000009 

0.53 

0.6017 

1.0451e-8 

1.861e-8 

0.56 

0.5821 

0.0000034 

0.000003 

1.04 

0.3129 

-0.000006 

0.000003 

-2.17 

0.0450 

Effect  Test 


Prediction  Profile 
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icreening  Fit 


Prediction  Profile 


•TSlIlllUl 


jcreening  Fit 


PAH  diss 


[Summary  of  Fit ) 

RSquare 

0.68452 

RSquare  Adj 

0.408476 

Root  Mean  Square  Error 

8375.271 

Mean  of  Response 

5796.961 

Observations  (or  Sum  Wgts) 

31 

^Analysis  of  Variance 
Source  DF  Si 

Model  14 

Error  16 

C  Total  30 


DF 

Sum  of  Squares 

Mean  Square 

14 

2435188423 

1.7394e8 

16 

1122322623 

70145164 

30 

3557511046 

■[Parameter  Estimates  ] 

Term 

Estimate 

Std  Error 

t  Ratio 

Prob>lt| 

Intercept 

6293.3322 

8243.35 

0.76 

0.4563 

Salinity 

-875.3955 

647.7484 

-1.35 

0.1954 

OR400 

0.0170025 

1.561499 

0.01 

0.9914 

Energy 

-56.02128 

297.912 

-0.19 

0.8532 

Loading 

34.982893 

83.96738 

0.42 

0.6825 

Salinity*Salinity 

15.107808 

16.97622 

0.89 

0.3867 

OR400*Salinity 

0.0250849 

0.014956 

1.68 

0.1129 

OR400*OR400 

0.0000002 

0.000081 

0.00 

0.9983 

Energy*Salinity 

5.0158159 

2.658816 

1.89 

0.0775 

Energy*OR400 

0.0109507 

0.005816 

1.88 

0.0780 

Energy*Energy 

0.1415734 

2.567392 

0.06 

0.9567 

Loading*Salinity 

-1.494825 

0.886272 

-1.69 

0.1111 

Loading*OR400 

-0.003833 

0.001939 

-1.98 

0.0655 

Loading*Energy 

-0.583095 

0.344661 

-1.69 

0.1101 

Loading*Loading 

0.118695 

0.285266 

0.42 

0.6829 

Effect  Test 


Prediction  Profile 


60904.54 


1681.491 


Salinity 


OR400 


Energy 


Loading 
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APPENDIX  IV 


Polycyclic  Aromatic  Hydrocarbon  Data 
Response  Surface  Experiments 


Client  ID  NA  NA  NA  Trial  1'  Trial  2'  Trials’  Trial  4'  Trials'  Trials’  Trial?’  Trials’  Trials’  Trial  10’  Trial  11’  Trial  12’ 

BattellelD  VU66PB  W44PB'  W53PB  UU24FIL  UU25FIL  UT54FIL  UT55FIL  UU26FIL  UU27FIL  UT56FIL  UT57FIL  UU11FIL  UU12FIL  UT68FIL  UT69FIL 

Sample  Volume  (L)  1.000  1.000  1.000  0.870  0.870  0.800  0.700  0.840  0.910  0.720  0.790  0.880  0.870  0.990  0.995 
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Energy  Density  Spectra  for  Flume 
Wind-Wave  Settings 
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FLUME  TESTS  1  THROUGH  5 


AVERAGE  PARTICLE  SIZE  RESULTS 


Battelle  Flume  Tank  Particle  Size  Results  ( 1998  ) 
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FLUME  TESTS  1  THROUGH  5 

TPH  RESULTS  IN  FILTERED  AND 
UNFILTER  WATERS 
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FLUME  TESTS  1  THROUGH  5 

PAH  RESULTS  IN  FILTERED  AND 
UNFILTER  WATERS 


FLUME  TEST  1  -FILTERED  &  UNFILTERED  WATERS 
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Surrogate  Recoveries 

Naphthalene-d8  4  &  64  13  &  17  & 

Phenanthrene-dIO  76  92  91  92 

Chtysene-d12  89  90  90  88 


FLUME  TEST  5  -  FILTERED  &  UNFLITERED  WATERS 


FLUME  TESTS  1  THROUGH  5 
MASS  BALANCE  RESULTS 


Battelle  Flume  Tests  1  through  5 


FINAL  MASS  BALANCE  TABLES 


Mass  Bitumen 
Added  (g) 

Floating  (g) 

Siink(g) 

Bathtub 

Ring(g) 

Suspended 

(g) 

Dissolved 

(g) 

TEST  1 

26280 

45 

17211 

5 

9009 

10 

TEST  2 

10366 

9453 

-27 

861 

74 

5 

TESTS 

10366 

8000 

2004 

29 

326 

6 

TEST  4 

10366 

7030 

1192 

1407 

730 

6 

TESTS 

10366 

3630 

-157 

1212 

5668 

12 

SpUl  Conditions 

Floating 

(%) 

Sunk  (%) 

Bathtub 
Ring  (%) 

Suspended 

(%) 

Dissolved 

(%) 

TEST  1 

low  energy, 
freshwater 

0.2 

65.5 

0.0 

34.3 

0.0 

TEST  2 

high  energy, 
brackish  water 

91.2 

-0.3 

8.3 

0.7 

0.0 

TESTS 

low  energy, 
brackish  water 

112 

19.3 

0.3 

3.1 

0.1 

TEST  4 

high  energy 
saltwater 

67.8 

11.5 

13.6 

7.0 

0.1 

TESTS 

high  energy, 
freshwater 

35.0 

B 

11.7 

54.7 

0.1 
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Battelle  -  latroscan  Data  Flume  Tests  1  through  5 


Appendix  VI 


TPH  and  PAH  Data  -  Battelie  Flume  Tests  1  through  5 


TESTl 

TPH  (TOTAL) 
(mg/L) 

GRAY 

TPH  (DISSOLVED) 
(mg/L) 

INST 

PAH  (TOTAL) 

(ug/L) 

PAH  (DISSOLVED) 

(iig/L) 

TIME 

(hrs) 

1 

5063.29 

2.30 

20303.6 

6.6 

2 

5328.95 

NA 

NA 

NA 

4 

5064.94 

NA 

NA 

NA 

6 

5197.37 

NA 

NA 

NA 

12 

5031.85 

NA 

NA 

NA 

24 

4200.00 

3.05 

16198.9 

5.8 

48 

3184.00 

NA 

NA 

NA 

72Repl(15cm) 

2081.08 

NA 

NA 

NA 

72Rep2(15cm) 

2062.50 

NA 

NA 

NA 

72Rep3(15cm) 

1776.47 

NA 

NA 

NA 

72Repl(35ciii) 

2187.50 

2.62 

NA 

5.0 

72Rep2(35cm) 

2211.76 

2.53 

NA 

NA 

72Rep3(35cm) 

2280.00 

2.44 

6285.3 

4.1 

72Repl(55cm) 

2617.45 

NA 

NA 

NA 

72Rep2(55cm) 

2535.21 

NA 

NA 

NA 

72Rep3(55cin) 

2518.52 

NA 

NA 

NA 

NA  -  not  analyzed 
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TPH  and  PAH  Data  -  Battelie  Flume  Tests  1  through  5 


TESTS 

TPH  (TOTAL) 
(mg/L) 

GRAY 

TPHpiSSOLVED) 

(mg/L) 

INST 

PAH  (TOTAL) 
(ug/L) 

PAH  (DISSOLVED) 
(ug/L) 

TIME 

(hrs) 

1 

1878.21 

1.37 

7090.1 

1.3 

2 

1549.67 

NA 

NA 

NA 

4 

741.75 

NA 

NA 

NA 

6 

630.00 

NA 

NA 

NA 

12 

361.14 

NA 

NA 

NA 

24 

157.30 

1.28 

491.9 

0.3 

48 

155.57 

NA 

NA 

NA 

72Repl(15cm) 

73.45 

NA 

NA 

NA 

72Rep2(15cm) 

72.50 

NA 

NA 

NA 

72Rep3(15cm) 

78.19 

NA 

NA 

NA 

72Repl(35cm) 

81.40 

1.55 

185.75 

0.2 

72Rep2(35cm) 

79.07 

1.68 

NA 

NA 

72Rep3(35cm) 

80.46 

1.64 

NA 

NA 

72Repl(55cm) 

77.53 

NA 

NA 

NA 

72Rep2(55cm) 

107.47 

NA 

NA 

NA 

72Rep3(55cm) 

83.72 

NA 

NA 

NA 

TPH  and  PAH  Data  -  Battelle  Flume  Tests  1  through  5 


TESTS 

TPH  (TOTAL) 
(mg/L) 

GRAY 

TPH(DISSOLVED) 

(mg/L) 

INST 

PAH  (TOTAL) 
(ug/L) 

PAH  (DISSOLVED) 
(ugL) 

TIME 

(hrs) 

1 

2381.72 

1.91 

6718.2 

1.2 

2 

2396.74 

NA 

NA 

NA 

4 

2144.44 

NA 

NA 

NA 

6 

2135.87 

NA 

NA 

NA 

12 

2134.41 

NA 

NA 

NA 

24 

2122.34 

2.12 

2483.8 

0.4 

48 

1552.69 

NA 

NA 

NA 

72Repl(15cm) 

1289.47 

NA 

NA 

NA 

72Rep2(15cm) 

1364.13 

NA 

NA 

NA 

72Rep3(15cm) 

1472.22 

NA 

NA 

NA 

72Repl(35cm) 

1366.48 

3.25 

273.87 

0.5 

72Rep2(35cm) 

1381.11 

3.19 

NA 

NA 

72Rep3(35cm) 

1385.96 

2.62 

NA 

NA 

72Repl(55cm) 

1467.74 

NA 

NA 

NA 

72Rep2(55cm) 

1478.02 

NA 

NA 

NA 

72Rep3(55cm) 

1548.94 

NA 

NA 

NA 

NA  -  not  analyzed 
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